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NAS3-27752, NASA AST AOI 14 
DESIGN AND TEST OF FAN/NACELLE MODELS 
QUIET HIGH-SPEED FAN 
DESIGN REPORT 


1.0 SUMMARY 

1.1 Purpose 

The primary objective of the Quiet High-Speed Fan (QHSF) program was to develop an 
advanced high-speed fan design that will achieve a 6 dB reduction in overall fan noise over a 
baseline configuration while maintaining similar performance. The program applies and 
validates acoustic, aerodynamic, aeroelastic, and mechanical design tools developed by NASA, 
US industry, and academia. The successful fan design will be used in an AlliedSignal Engines 
(AE) advanced regional engine to be marketed in the year 2000 and beyond. This technology is 
needed to maintain US industry leadership in the regional turbofan engine market. 

1.2 Conclusions 

The AE/QHSF was designed to match the baseline rotor design point performance and yet 
achieve a 6 dB reduction. Innovative techniques in coupling the aerodynamic and mechanical 
analyses, along with design of experiments (DOE’s) involving multiple engineering disciplines 
accomplished the effort. The result was a fan stage design with a forward swept rotor with 50 
degrees of leading edge sweep that matched the baseline aerodynamic performance and achieved 
all mechanical and aeroelastic design criteria; and, a stator lean tailored for minimum rotor-stator 
acoustic interaction. 

1.3 Design Procedure 

The design technique used for the QHSF involved early design input from the Aerodynamic, 
Mechanical, Acoustic, and Aeroelastic disciplines through the use of design of experiments. 
Additional optimization procedures were used to enhance and verify the DOE analyses and 
results. Four DOE’s were completed for the rotor and two for the stator which focused on 
achieving the design goals without violating established design criteria. 
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1.4 Design Summary/Perspective 


In order to obtain the design goal of a 6 dB reduction, the design concentrated on three aspects. 
Reduction of the inlet relative normal Mach number, position of the shock within the passage, 
and minimizing rotor-stator interaction effects. Leading edge effective sweep was used to reduce 
the inlet relative normal Mach number, while mean-camber line blade angle distributions were 
used to correctly position the shock in the blade-to-blade passage. Although some design 
constraints were imposed on the stator design, rotor-stator interaction effects were minimized by 
controlling the rotor wake impingement on the stator leading edge. 

The AE/QHSF is a forward swept rotor designed to match the design point performance of the 
baseline. It incorporates relatively large amounts of leading edge effective sweep, and yet with 
the use of statistical analyses and optimization procedures has reduced stress and weight while 
maintaining adequate aeroelastic and critical vibration margins relative to the baseline. Table 
1.4-1 shows a summary of the AE/QHSF design features at the aerodynamic design point. 
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TABLE 1.4-1 AE/QHSF SUMMARY 
AERODYNAMIC DESIGN POINT 


Performance Parameters: 


Inlet corrected flow 98.65 lbm/s 

Bypass ratio 3.78 


DAWES Results in AXCAPS Test 

QHSF Baseline Baseline 

Rotor Stage Rotor Stage Stage 


Bypass pressure ratio 1.839 1.813 1.839 1.811 1.844 

Core pressure ratio 1.845 1.817 1.777 1.741 1.718 

Bypass efficiency (ad) 0.883 0.860 0.888 0.865 0.889 

Core efficiency (ad) 0.960 0.934 0.961 0.923 0.918 

O/A pressure ratio 1.813 1.796 1.818 

O/A efficiency (ad) 0.874 0.875 0.895 

Rotor Parameters: 

Number of blades 22 

Corrected RPM 15358 

Inlet tip diameter 22.0 inches 

Corrected tip speed 1474 ft/s 

Hub - tip ratio 0.35 

Aspect ratio 1.35 

Inlet specific flow 42.6 lbm/s/ft 2 

Hub and tip slope 28°, -7° 

Hub work coefficient 2.5 

Inlet tip rel Mach no 1.43 

Stator Parameters: 

Number of vanes 52 

Aspect ratio 2.42 

Exit Mach no. (avg.) 0.52 

Exit air angle (avg.) 0° 
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2.0 


INTRODUCTION AND BACKGROUND 


2.1 AST Area of Interest 14 

Area of Interest 14 specifically addresses two level-three milestones to support "Model Tests for 
Code Validation/Concept Evaluation" in the Noise Reduction element of the NASA Advanced 
Subsonic Technology (AST) Program: 

• Select second- generation, low-noise concepts for model tests (4Q, FY 1997) 

• Complete second-generation model tests for low noise designs (4Q, FY 1999) 

A 2-year effort was originally proposed for this effort. During the first year, a calibration of the 
acoustic and aerodynamic prediction methods was performed and a baseline fan definition was 
established and evaluated. The second year activities were to include evaluation of several 
candidate noise reduction concepts using combinations of rotor and stator lean and sweep. The 
program scope has since been expanded. The program roadmap is illustrated in Figure 2.1-1. 
Three fan designs are used in the program. The existing NASA-Fewis QF-12 fan and an AE 
TFE derivative fan serve as baseline configurations. Data from these fans will aid in design tool 
validation. A new advanced fan was designed by AE. The program will lead to a model scale 
test using the Ultra High Bypass (UHB) Propulsion Simulator in the NASA-Lewis 9 x 15-foot 
wind tunnel. 


Fan 


QF-12 


Baseline 


AlliedSignal 

Advanced 

Design 


1995 

Calibration of 
Dawes CFD 
Code 


Calibration of 
Dawes CFD & 
UTRC/ 
Eversman 
Acoustic Code 


Optimum Stator 
Aerodynamic 
Design 


1996 


1997 


Fabricate Fan 


Optimum Rotor 
Design / Final 
Fan Aero 
Design 


Mechanical and 
Aeroelastic 
Analysis 


1998 


Baseline for 
Rig Test 


Fabricate Fan / 
Rig Test 


1999 


Data Reduction 
/ Data Analysis 

/ 

Reporting 

Results 


Figure 2.1-1. The Quiet High-Speed Fan Program Roadmap Targets a Wind 

Tunnel Test in October of 1998. 
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2.2 Quiet High Speed Fan Program 


The QHSF program is a cooperative effort between AE and the NASA Lewis Research Center. 
Both NASA and AE will design an advanced high-speed fan that will be tested on the Universal 
Propulsion Simulator in the NASA Lewis 9 x 15 foot wind tunnel, currently scheduled for the 
third quarter of 1998. An AE scaled TFE derivative fan design will be used as a baseline. A 
nacelle model will be provided that will be characteristic of a typical, modern regional aircraft 
nacelle and meet all of the program test objectives. 

The quiet high-speed fan is an advanced single-stage fan designed for a 5K to 20K pound thrust 
turbofan regional airline application. Advanced aerodynamic, mechanical, aeroelastic, and 
computational fluid dynamic (CFD) tools will be used to meet aggressive performance goals 
while achieving at least a 6 dB reduction in fan noise at a critical takeoff noise condition. Three 
fans will be built for evaluation at NASA; a baseline TFE derivative, an advanced NASA design, 
and an advanced AE design. A 22-in. diameter fan nacelle for the UHB Propulsion Simulator in 
the NASA-Lewis 9x15 foot wind tunnel will integrate NASA and AE aerodynamic designs and 
instrumentation to fulfill all aerodynamic, mechanical, distortion sensitivity, and acoustic test 
needs. 


2.3 Baseline Fan 

The baseline consists of a damperless, low-aspect- ratio, moderately aft swept rotor and full-span 
aft swept composite stator vanes. The geared fan configuration allows the fan to run at a tip 
speed optimized for best performance, stall margin, and noise. 

The baseline fan was the subject of considerable acoustic evaluation early in the design phase. 
Blade and vane counts, rotor/stator spacing, and stator vane sweep were selected to minimize the 
noise signature within the fan design constraints. Fan rig spinning-mode measurements were 
made to verify the design constraints and provide acoustic treatment design criteria. The acoustic 
effort was concluded with a full-scale engine acoustic test which verified that the acoustic design 
goals, including the acoustic treatment, were satisfied. 

Installed performance of the baseline fan was a critical issue. Measuring fan component 
performance accurately during flight was required. The engine nacelle and inlet section of the 
flight test engine were custom tailored and instrumented for flow measurement. The engine front 
frame was also instrumented so fan performance could be measured without impacting fan or 
engine performance. The fan performance measured in flight at appropriate Reynolds numbers 
during the extensive flight test program agreed very well with rig and engine data acquired on the 
ground using conventional fluid metering techniques. The experience gained on these programs 
will be used on the Quiet High-Speed Fan program. 
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2.4 Historical Background 


The presence of a shock at the inlet of the fan rotor in a turbofan engine can result in acoustic 
phenomena that represent substantial noise sources. Multiple Pure Tone (MPT) noise, for 
example, results when the pressure disturbances from the inlet shock move upstream out of the 
rotor blade passage. One approach to reducing these shock-related noise sources is to eliminate 
the formation of the inlet shock in the fan by tailoring the rotor blade shape. The introduction of 
sweep in the fan rotor blade can reduce the relative velocity component normal to the blade to 
subsonic values, much as a swept wing on an aircraft can produce subsonic velocities normal to 
the wing leading edge, even when the resultant velocity is supersonic. 

This noise minimization technique was applied by AlliedSignal Engines (formerly AVCO 
Lycoming) and Bolt, Beranek and Newman, Inc. (BBN) to the design of the QF-12 quiet high- 
speed fan, as part of a NASA-sponsored program performed between 1974 and 1977 (Reference 
1). The fan rotor featured a compound forward- and-aft sweep to eliminate the leading edge 
shock. However, although shock-induced MPT noise was reduced, the aerodynamic 
performance of the fan did not meet design goals. The reasons for the performance deficiencies 
were not determined, but the source was found to be localized in the rotor. Figure 2.4-1 shows a 
summary of the results of the acoustic evaluation of the QF-12. 

The ability to accurately predict, during the fan rotor design, the speed range over which the MPT 
noise will occur can provide a valuable tool for acoustically tailoring the design of the rotor to 
minimize the effect of this noise source. 



Figure 2.4-1. The QF-12 Fan Demonstrated the Feasibility of Using Blade 

Leading Edge Sweep to Reduce Fan Noise. 
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3.0 


DESIGN OBJECTIVES / CRITERIA 


3.1 Overall 

The AE/QHSF stage design included design objectives from the acoustic, aerodynamic, 
mechanical, and aeroelastic disciplines. Acoustic objectives included a 6 dB reduction in 
effective perceived noise (EPNdB) at takeoff relative to the baseline. This was to be 
accomplished without the known benefits of vane count reduction and increased axial spacing 
(hub and tip) between rotor TE and vane LE. The intent of the design was focused on obtaining 
noise reductions through the use of more unconventional design techniques such as: controlling 
the rotor wake impinging on the stator LE, reducing rotor normal inlet relative Mach number, 
and controlling shocks within blade passage. Aerodynamic objectives were to obtain baseline 
performance (or better) with a forward swept rotor and acoustically matched vanes. Mechanical 
and aeroelastic objectives were to obtain AE and NASA criteria for stress and burst margins, 
proper positioning of mode / harmonic crossings to minimize high vibratory strains, and provide 
sufficient flutter margin throughout the operating region. 


3.2 Specific 

Specific aerodynamic design criteria is shown in Table 3.2-1. Mechanical criteria is shown in 
Table 3.2-2 and Figure 3.2-1. 


Table 3.2-1. Specific Aerodynamic Design Criteria 


Fan Aero Desian Point (100% Nc) 



22” DIA 

Wcorr, Ibm/s 

98.9 

Wc/A, lbm/s/ft A 2 

42.7 

Utcorr, ft/s 

1474 

Bypass ratio 

3.8 

P/P, overall 

1.82 

Eff ad, overall 

> .895 

Stall marciin (Nio=C) 

15% 

Hub/Tip ratio 

.35 

Rotor blade count 

22 

Stator vane count 

52 
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Table 3.2-2. Specific Mechanical Design Criteria 


• Airfoil and Attachment LCF Life > 1 0 A 5 cycles 

- Peak Stress in Attachment and Airfoil 

• HCF Life > 1 0 A 7 Cycles 

- Frequency, Excitation (Inlet Distortion Harmonics, Speed 
(©Resonance) 

• Bird-Ingestion Capability 

- Bird Weight, LE Thickness 

• Flutter Margin Exceeds Baseline 

- Reduced Frequency, Twist/Flex, and Incidence 

• Blade Weight Equal to Or Less Than Baseline 


• High Speed Resonances Must Be 
Avoided 


- Rotor 

• Mode 1 / 2E Crossing < 50 
percent speed 

• Mode 2 / 4E Crossing < 70 
percent speed 

10% Frequency Margin at 
Maximum Engine Operating Speed 
Modified Goodman Approach used 
for HCF Calculation 



Figure 3.2-1. Mechanical Criteria for High-Cycle Fatigue 
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4.0 PRELIMINARY DESIGN 


4.1 Cycle Point Selection 

An engine cycle design point was originally established for the AE/QHSF to correspond to a 
regional type application for the AE Small Engine Technology (SET) engine concept. However, 
during the initial study efforts, it became apparent that although the cycle point originally 
selected would benefit the AST program in terms of a regional application; it would make 
comparisons relative to the baseline (design and test) much more difficult. Therefore, the 
baseline design cycle was selected which would provide the best and cleanest ‘back-to-back’ test 
to verify Area of Interest 14 technology goals. 

4.2 Rotor 

4.2.1 Blade Harmonic Selection 

Analyses were completed to determine whether the rotor should be a “1/rev” design (similar to 
the baseline) or a “2/rev”. The 1/rev design offers the advantage of less weight but, without a 
midspan damper, increases the susceptibility to flutter. An optimization process using ANSYS 
was completed for both designs in order to minimize weight while maintaining all other design 
criteria. 

A parametric model of the airfoil was defined to allow aeromechanical optimization via 
ANSYS. The value of maximum airfoil thickness (Tmax/C) and its chordwise position were 
varied at several spanwise locations, including the hub and tip. A cubic spline fit through these 
control points was used to define these distributions at 5% span increments. 

X-cg and Y-cg distributions were defined along 4 th -order polynomial functions fixed at the airfoil 
hub. The X-cg at the tip was fixed at 1.07 inches forward. The coefficients of the polynomials 
were used as design variables. X-cg and Y-cg were defined at 5% span increments. 

The blade geometric parameters were used to create an AXCAPS (AE’s radial equilibrium 
streamline flow solver) input file, with the resulting airfoil geometry automatically meshed and 
analyzed. These geometric parameters were defined as independent design variables for the 
optimization procedure. A single objective function is available with ANSYS, thus additional 
“objectives” were forced with constraints. 

The first optimization was for a blade with a 1/rev crossing (1E/4E). The initial analyses were 
conducted with the Tmax/C and position of Tmax design variables defined at 4 equally-spaced, 
spanwise sections. The resulting optimized airfoil was approximately 23 percent heavier than the 
baseline design. The spanwise resolution was subsequently increased to 6 sections. A cubic 
spline fit through these control points was used to define these distributions at 5% span 
increments. Using the new spanwise resolution, the optimized blade was only 1 percent heavier 
than the baseline blade. 
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Table 4.2. 1-1. AE/QHSF Preliminary 1/rev Design 


Parameter 

QHSF 1/Rev Airfoil 
(Fixed-root; no root fillet) 

Relative to Baseline 

Airfoil Weight 

1.27 lb 

+ 1.2 % 

Max Principal Stress 

82 ksi 

-10.8 % 

Max Deflection 

0.31 in 

-19.5 % 

Mode 1/2E crossing 

48 percent speed 

-28 % 

Flutter Margin 

-0.44 

+10% 

Tip Tmax/C 

0.02 

Same 

Initial optimization analyses for the 2/rev crossing (2E/4E) were also conducted with Tmax/C 
and position of Tmax design variables defined at 4 equally-spaced, spanwise sections. The best, 
near-feasible solution obtained was 31 percent heavier than baseline blade. Based on the success 
of the 1E/4E study, the spanwise Tmax/C and position of Tmax design variables were defined up 
to 11 sections. A cubic spline fit through these control points was used to define these 
distributions at 5% span increments. The total number of design variables (DY’s) was 29 for the 

11 point case even though ANSYS recommends no more than 

20 DV’s. 

The weight of the 2E/4E design was relatively unaffected by 

an increase in spanwise sections. 

The best design was slightly low 

on mode 1 frequency margin, 14 percent vs. the desired 15 

percent, but met all other mechanical constraints. The blade was still 22 percent heavier than the 
baseline blade; and, was obtained only after reducing Tmax/C at the tip to 0.01, well below the 

design criteria of 0.02. 



Table 4.2. 1-2. AE/QHSF Preliminary 2/rev Design 

Parameter 

QHSF 2/Rev Airfoil 
(Fixed-root; no root fillet) 

Relative to Baseline 

Airfoil Weight 

1.53 lb 

+22% 

Max Stress 

67 ksi 

-27% 

Max Deflection 

0.39 in 

0 

Mode 2/4E freq margin 

14 percent (require 15%) 


Flutter Margin 

-0.32 

+18% 

Tip Tmax/C 

0.01 

-50% 
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4.2.2 DAWES Model 


The DAWES program is a finite- volume time-marching solver for the 3-D thin-layer Navier- 
Stokes equations. The version of the program used by AE for the QHSF was the same version 
previously used on the baseline rotor design. To maintain close compatibility, the same grid and 
flow size were used for the two designs. All DAWES analyses, rotor and stator, were completed 
at the baseline engine flow size. 

The blade and flowpath geometry descriptions for the fan rotor were obtained from AE’s 
axisymmetric streamline curvature program, AXCAPS. This program provides a discrete-point 
definition of the blade geometry along specified stream surfaces intersected with the blade 
surface. Endwall definition is also provided in terms of radial and axial coordinates at discrete 
points. The AXCAPS model of the flowpath is shown in Figure 4.2.2- 1. 



Axial Length, in. 

Figure 4.2.2-1. AXCAPS Flowpath Model Used 


The geometry represented in the AXCAPS model was the deflected aerodynamic blade shape at 
100 percent wheel speed. Hub fillets were not modeled in DAWES; and, the blade tip was 
modeled as tapered, rather than squared-off, which is required by a gridding limitation within the 
program. 
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The computational domain consisted of a single rotor blade passage, with the suction surface of 
one blade forming one pitchwise boundary, and the pressure surface of the adjacent blade 
forming the other pitchwise boundary. The inlet to the computational domain was positioned 
upstream (approximately one tip chord) to reduce impacting the development of the upstream 
shock structure. The exit of the computational domain was located the same, relative to the rotor 
TE, as the baseline. Because of the rotor forward sweep, the exit plane did not include the outer 
40% span of the stator leading edge. Analyses were completed late in the design with an 
extended grid which showed there was no impact on the previously completed rotor wake/stator 
interaction study. 


The computational grid for the DAWES program is a structured single-block skewed H-grid. 
Typical sections of the computational grid are presented in Figure 4.2.2-2. The grid size 
employed for the final set of rotor analyses was specified as: 

Pitchwise: 25 Nodes 

Spanwise: 71 Nodes (with 8 cells in the clearance gap) 

Streamwise: 131 Nodes 
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R vs. R*phi 
mid-passage 


Figure 4.2.2-2. Rotor DAWES Computational Grid 
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4.2.3 Design of Experiment 1 
4.2.3. 1 Purpose 

The first statistical design of experiment (DOE) was created to determine the extent of forward 
sweep required to obtain the 6 dB noise reduction with marginal violation of the other 
performance goals. (It was believed that the performance goals would be achieved in subsequent 
DOE’s planned.) The DOE involved analyses from the Aerodynamic, Mechanical, Acoustic, and 
Aeroelastic disciplines. 


4.2.3.2 Description 

Because of the perceived non-linearity of the axial and tangential leans, the DOE was created to 
consist of 4 factors (Xcg-tip, Ycg-tip, Xcg-midspan, Ycg-midspan) with 3 levels each. The 9 
cases are described in Table 4.2. 3-1, with the CG orientation shown in Figure 4.2.3-1. Figures 
4.2.3-2 and 4.2. 3-3 show the CG distribrutions of the nine cases. The blade geometry and the 
mean flow properties were determined using AXCAPS. The intent was to analyze the cases at 
design speed (100% Nc) and 80% Nc using DENTON (Denton’s 3D inviscid Euler code) to 
obtain the aerodynamic and acoustic quality characteristics. It was believed, based on past 
experience, that DENTON could be used to provide accurate solutions for flow, and shock 
locations; and, because of the fast turnaround time would be advantages for the DOE. For 
solution verification, several cases were analyzed using DAWES which showed some differences 
in flow and passage shock locations. The differences were believed to be caused by using 
convergence factors (time step, pitch and axial smoothing, damping, and inlet smoothing) in 
DENTON which deviated from experience but were required for case convergence. DAWES 
was therefore used to analyze the cases; however, because of the CPU time required, only the 
design speed was analyzed for the DOE. 

Table 4.2.3-1. DOE 1 Forward Sweep Analysis X-Factors 


Case No. 

X-cg tip 

Y-cg tip 

X-cg mid 

Y-cg mid 

1 

-1.50 

-1.0 

-0.25 

-0.50 

2 

-1.50 

0 

0 

0 

3 

-1.50 

+1.0 

+0.60 

+0.50 

4 

-0.75 

-1.0 

0 

+0.50 

5 

-0.75 

0 

+0.60 

-0.50 

6 

- 0.75 

+1.0 

-0.25 

0 

7 

0 

-1.0 

+0.60 

0 

8 

0 

0 

-0.25 

+0.50 

9 

0 

+1.0 

0 

-0.50 


Dimensions are in inches and relative to Baseline engine size (30.7 in Dia.). 
Negative X-cg is axial forward direction, positive Y-cg is in direction of 
rotation 
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Figure 4.2.3-1. Blade CG Orientation Definition. 


Baseline 



Figure 4.2.3-2. Xcg distributions used in DOE 1. 
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Baseline 



Figure 4.2.3-3. Ycg distributions used in DOE 1. 


4.2.3.3 Analysis 

Response variables (Y-factors) were compiled based on analyses from the four engineering 
disciplines. These factors became somewhat dynamic in that some were deleted, while others 
were added or calculated slightly differently. The final list of Y-factors with a brief description is 
shown in Table 4.2. 3-2 with results tabulated in Table 4. 2.3-3. Vector solutions for each case 
were created at 55.9% corrected fan speed using an off-design model from the design point. This 
speed corresponds to an approach condition where interaction noise is likely to be the dominant 
fan noise source (buzzsaw noise cutoff). Predictions of rotor/stator interaction noise were 
performed using the NASA code V072. The baseline stator was used for all rotor design studies. 
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Table 4.2.3-2. DOE 1 Forward Sweep Analysis Y-Factors 



Discipline 

Y-factor 

Description 

Analvsis Tool 

% Nlc 

1 

Acoustic 

dB 

Interaction noise sound power level 

V072 

55.9 

2 

Acoustic 

eff swp le 

Rotor LE effective sweep 

AXCAPS 

100.0 

3 

Acoustic 

%c 50 

Passage shock location at 50% span 

DAWES 

100.0 

4 

Acoustic 

%c 60 

Passage shock location at 60% span 

DAWES 

100.0 

5 

Acoustic 

%c 70 

Passage shock location at 70% span 

DAWES 

100.0 

6 

Acoustic 

%c 80 

Passage shock location at 80% span 

DAWES 

100.0 

7 

Acoustic 

%c 90 

Passage shock location at 90% span 

DAWES 

100.0 

8 

Acoustic 

%c 95 

Passage shock location at 95% span 

DAWES 

100.0 

9 

Acoustic 

delz 50 

Axial distance RTE - VLE midspan 

AXCAPS 


10 

Acoustic 

delz 100 

Axial distance RTE - VLE tip 

AXCAPS 


11 

Aero 

Wc 

Inlet corrected flow (top of choke) 

DAWES 

100.0 

12 

Aero 

PR 

Rotor pressure ratio (top of choke) 

DAWES 

100.0 

13 

Aero 

Eff 

Rotor efficiency (top of choke) 

DAWES 

100.0 

14 

Aero 

Mpeak 50 

Suction surface peak Mach 50% span 

DAWES 

100.0 

15 

Aero 

Mpeak 60 

Suction surface peak Mach 60% span 

DAWES 

100.0 

16 

Aero 

Mpeak 70 

Suction surface peak Mach 70% span 

DAWES 

100.0 

17 

Aero 

Mpeak 80 

Suction surface peak Mach 80% span 

DAWES 

100.0 

18 

Aero 

Mpeak 90 

Suction surface peak Mach 90% span 

DAWES 

100.0 

19 

Aero 

Mpeak 95 

Suction surface peak Mach 95% span 

DAWES 

100.0 

20 

Aero 

ws 50 

Wennerstrom shock loss 50% span 

AXCAPS 

100.0 

21 

Aero 

ws 75 

Wennerstrom shock loss 75% span 

AXCAPS 

100.0 

22 

Aero 

ws 100 

Wennerstrom shock loss 100% span 

AXCAPS 

100.0 

23 

Aero 

spds 

Speedline slope out of choke 

DAWES 

100.0 

24 

Aeroelastic 

flut 1 

Empirical flutter parameter mode 1 

ANSYS 

100.0 

25 

Aeroelastic 

flut 2 

Empirical flutter parameter mode 2 

ANSYS 

100.0 

26 

Aeroelastic 

twf 1 

Twist/flex ratio mode 1 

ANSYS 

100.0 

27 

Mechanical 

ssavgl 

Avg suction surface stress 4 - 20% 

ANSYS 

100.0 

28 

Mechanical 

ssavg2 

Avg suction surface stress 21 - 40% 

ANSYS 

100.0 

29 

Mechanical 

ssavg3 

Avg suction surface stress 41 - 60% 

ANSYS 

100.0 

30 

Mechanical 

ssavg4 

Avg suction surface stress 61 - 80% 

ANSYS 

100.0 

31 

Mechanical 

ssavg5 

Avg suction surface stress 81 - 100% 

ANSYS 

100.0 

32 

Mechanical 

psavgl 

Avg pressure surface stress 4 - 20% 

ANSYS 

100.0 

33 

Mechanical 

psavg2 

Avg pressure surface stress 21 - 40% 

ANSYS 

100.0 

34 

Mechanical 

psavg3 

Avg pressure surface stress 41 - 60% 

ANSYS 

100.0 

35 

Mechanical 

psavg4 

Avg pressure surface stress 61 - 80% 

ANSYS 

100.0 

36 

Mechanical 

psavg5 

Avg pressure surface stress 81 - 100% 

ANSYS 

100.0 

37 

Mechanical 

umax tip 

Max tip deflection 

ANSYS 

100.0 

38 

Mechanical 

freq 1 

Frequency margin mode 1 

ANSYS 

100.0 

39 

Mechanical 

freq 2 

Frequency margin mode 2 

ANSYS 

100.0 

40 

Mechanical 

freq 3 

Frequency margin mode 3 

ANSYS 

100.0 

41 

Mechanical 

fee 1 

Placement of 2/rev crossing mode 1 

ANSYS 

100.0 

42 

Mechanical 

fee 2 

Placement of 4/rev crossing mode 2 

ANSYS 

100.0 
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The process flowchart used for DOE 1 is shown in Figure 4.2.3-4. The flowchart shows how the 
traditional “experiment” analysis is coupled with the ANSYS optimization technique. 



Figure 4.2.3-4. The process flow diagram for the QHSF design shows the 

combined use of design of experiments and ANSYS parametric 
optimization. 


After all the response variables were collected, the data was input into the statistical program 
JMP. A second order polynomial was used to fit the data, and coefficients were generated for 
each Y-factor. Because of some deficiencies in the JMP program, the coefficients were input 
into EXCEL and used with the SOLVER function. SOLVER works by maximizing, minimizing, 
or iterating to a specified value for one Y-factor and can extrapolate outside the X-factor 
envelope. To obtain a feasible solution, limits were set for the other Y-factors and weighting was 
somewhat controlled by the range of these limits. Various feasible solutions were obtained by 
optimizing different Y-factors. Several new rotor geometry cases (go-forward cases) were 
created based on these results and analyzed. Go-forward cases created using EXCEL Solver and 
the statistical program JMP were designated jmp(number), while cases created from the ANSYS 
optimization were designated opt(letter). 
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4.2.3.4 Results 


Acoustic response variables, overall power level at 55.9% Nlc relative to the baseline rotor and 
blade effective sweep at the leading edge and at the suction surface impingement point of the 
adjacent blade, are shown in Figure 4.2. 3-5. Although the figure shows only one case (case #4) 
to have lower noise than the baseline, the purpose of the DOE was to acquire sensitivities of the 
design variables and response variables. The sensitivities were the objective of the experiment 
such that go-forward cases could be created. The figure also shows two important aspects of 
effective sweep: the large effect of tangential lean, and the importance of the suction surface 
impingement point. Figure 4.2.3-2 showed that only one case (case #8) had a small amount of aft 
sweep; however Figure 4.2. 3-5 shows that cases with large amounts of tangential lean (suction 
side down) can cause the effective sweep to be in the opposite direction thereby (potentially) 
mitigating benefits of the forward swept blade. 

The blade effective sweep was calculated based on Art Wennerstroms’ oblique shock angle 
(OB A) calculation at the leading edge and the shock impingement on the adjacent airfoil’s 
suction surface. In AE’s coordinate system, oblique shock angles greater than 90 degrees 
represent forward sweep, and are shown in the figure as a negative value (relative to 90 degrees). 
The desired amount of effective sweep was based on the amount required to reduce the normal 
component of the inlet relative Mach number at 90% Nlc (maximum obtainable fan speed up to 
4000 feet elevation) to a Mach number of 1.0. Effective sweep described here is the delta 
between the angle of the shock surface relative to the engine center line and 90 degrees. The 
shock surface at the leading edge is created as the shock travels from the leading edge of one 
blade and impinges on the adjacent blade’s suction surface at some distance aft of the leading 
edge. Therefore, the shock loss (at a particular radius) is a function of the oblique shock angle at 
the leading edge and the impingement point on the suction surface. Because the shock surface 
impingement point on the adjacent blade is aft of the blade leading edge, the oblique shock angle 
at the impingement point will always be less than the oblique shock angle at the leading edge. 
Therefore, for aft swept and forward swept blades with the same absolute value of leading edge 
effective sweep, the aft swept blade will have lower shock losses. 

Aerodynamic response variables, inlet corrected flow and rotor efficiency, are shown in Figure 
4.2.3-6. Several empirical aeroelastic response variables for Mode 1, twist/flex ratio and 
empirical flutter parameter (EFP), are shown in Figure 4.2.3-7. The EFP is the ratio of the actual 
reduced frequency and the desired reduced frequency for a twist/flex ratio, based on AE’s 
empirical data base. Figure 4.2.3-8 shows two, of many, mechanical response variables. The 
two variables shown are the maximum tip deflection and the average bending stress between 41 
and 60% radial span. The tip deflection is the vector summation of the radial, axial, and 
tangential components and is relative to the baseline engine size. The goal of +/- 0.4 inches 
shown is also relative to the baseline engine size. The dominant component of the deflection 
vector is axial, while the radial component is much less (approximately 0.060 inches). Positive 
values of tip deflection indicate bending toward the suction side. The bending stress goal of less 
than +/- 40 ksi on the suction and pressure surfaces was used to keep the maximum principal 
stresses on the blade below 85 ksi. 
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Overall Power Level 
Relative to Baseline 

Blade Effective Sweep 
(degrees) 


2.5 


2 

1.5 

1 

0.5 



1 23456789 Baseline 

Case No. 


Figure 4.2.3-5. Result of the Acoustic Assessment of the Configurations for 
Rotor DOE 1. 
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Figure 4.2.3-6. 
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Figure 4.2.3-7. Result of the Aeroelastic Assessment of the Configurations for 
Rotor DOE 1. 
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Figure 4.2.3-8. Result of the Mechanical Assessment of the Configurations for Rotor 
DOE 1. 
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Table 4.2.3-3. Tabulated Results of DOE 1 Forward Sweep Analysis 



Config # 

i 

2 

3 

4 

5 

6 

7 

8 

9 

1 

dB 

132.4 

132.8 

133.0 

131.4 

133.8 

132.1 

133.1 

132.8 

132.1 

2 

eff swp le 

-6.27 

-14.70 

-23.53 

22.61 

-25.50 

-23.57 

8.73 

10.28 

-28.40 

3 

%c 50 

28.80 

31.26 

33.66 

35.85 

31.00 

26.49 

35.79 

30.91 

21.83 

4 

%c 60 

33.61 

36.29 

39.09 

43.65 

36.19 

33.56 

46.53 

38.36 

28.30 

5 

%c 70 

41.13 

47.06 

47.42 

48.86 

44.33 

43.87 

52.07 

46.42 

38.11 

6 

%c 80 

46.65 

55.76 

59.04 

51.82 

53.06 

52.46 

60.46 

52.40 

46.26 

7 

%c 90 

55.96 

66.79 

69.56 

55.31 

66.93 

68.87 

68.49 

55.96 

65.99 

8 

%c 95 

58.55 

73.37 

75.69 

55.48 

73.39 

73.26 

70.63 

58.83 

70.62 

9 

delz 50 

3.827 

3.601 

3.024 

3.526 

3.004 

3.757 

2.944 

3.767 

3.516 

10 

delz 100 

6.862 

6.784 

6.737 

6.228 

5.950 

6.036 

5.298 

5.465 

5.235 

11 

Wc 

191.01 

191.69 

189.59 

190.39 

189.94 

188.99 

191.02 

191.16 

187.67 

12 

PR 

1.899 

1.901 

1.894 

1.811 

1.898 

1.890 

1.852 

1.854 

1.891 

13 

Eff 

.9038 

.9082 

.8999 

.8945 

.9054 

.9010 

.9052 

.9031 

.8942 

14 

Mpeak 50 

1.429 

1.464 

1.470 

1.470 

1.453 

1.404 

1.480 

1.422 

1.423 

15 

Mpeak 60 

1.452 

1.480 

1.491 

1.496 

1.483 

1.432 

1.504 

1.441 

1.440 

16 

Mpeak 70 

1.454 

1.464 

1.475 

1.503 

1.468 

1.432 

1.545 

1.444 

1.448 

17 

Mpeak 80 

1.445 

1.460 

1.475 

1.514 

1.459 

1.426 

1.550 

1.454 

1.429 

18 

Mpeak 90 

1.467 

1.527 

1.535 

1.510 

1.554 

1.545 

1.574 

1.517 

1.503 

19 

Mpeak 95 

1.493 

1.534 

1.532 

1.465 

1.586 

1.631 

1.571 

1.488 

1.613 

20 

ws 50 

.023 

.028 

.034 

.016 

.029 

.024 

.019 

.019 

.026 

21 

ws 75 

.055 

.046 

.031 

.045 

.043 

.040 

.053 

.046 

.039 

22 

ws 100 

.081 

.076 

.055 

.040 

.048 

.045 

.077 

.059 

.029 

23 

spds 

-9.5e-4 

2.79e-3 

3.90e-3 

-2.84e-3 

3.01e-3 

9.34e-3 

-1.09e-3 

7.35e-4 

-9.5e-4 

24 

flut 1 

-0.563 

-0.546 

-0.568 

-0.627 

-0.501 

-0.473 

-0.594 

-0.516 

-0.448 

25 

flut 2 

-0.105 

-0.008 

-0.129 

-0.181 

-0.005 

-0.005 

-0.008 

0.319 

-0.200 

26 

twf 1 

0.581 

0.562 

0.581 

0.851 

0.533 

0.484 

0.696 

0.561 

0.468 

27 

ssavgl 

152.5 

71.9 

-14.5 

92.2 

21.4 

-1.3 

48.1 

21.2 

-50.3 

28 

ssavg2 

221.0 

164.1 

110.5 

232.5 

49.6 

-25.5 

120.6 

54.8 

-152.9 

29 

ssavg3 

235.1 

223.6 

250.7 

296.9 

115.8 

-21.4 

169.4 

54.0 

-205.0 

30 

ssavg4 

173.5 

184.0 

241.9 

200.6 

113.6 

-12.1 

149.1 

28.4 

-138.9 

31 

ssavg5 

39.5 

45.1 

63.4 

48.3 

31.8 

-3.3 

38.0 

4.0 

-35.7 

32 

psavgl 

-110.6 

-14.8 

92.9 

-45.5 

30.6 

69.8 

3.4 

44.1 

108.7 

33 

psavg2 

-191.9 

-126.8 

-66.3 

-207.8 

-6.5 

82.0 

- 87.6 

-9.8 

213.9 

34 

psavg3 

-211.6 

-195.9 

-221.4 

-283.4 

-81.2 

65.3 

-150.1 

-22.5 

256.3 

35 

psavg4 

-152.1 

-160.0 

-215.2 

-192.9 

-88.6 

38.2 

-133.8 

-6.9 

165.3 

36 

psavg5 

-32.4 

-38.1 

-57.0 

-45.2 

-24.5 

11.8 

-31.7 

3.6 

44.1 

37 

umax tip 

-4.83 

-3.64 

-3.05 

-4.35 

-1.48 

1.22 

-2.50 

-0.35 

4.78 

38 

freq 1 

1.478 

1.491 

1.461 

1.478 

1.559 

1.504 

1.606 

1.586 

1.527 

39 

freq 2 

3.541 

3.630 

3.448 

3.242 

3.775 

3.751 

3.640 

3.827 

3.166 

40 

freq 3 

0.692 

0.533 

0.802 

0.785 

0.423 

0.589 

0.498 

0.304 

0.958 

41 

fee 1 

0.420 

0.457 

0.436 

0.423 

0.556 

0.518 

0.579 

0.610 

0.407 

42 

fee 2 

0.724 

0.784 

0.673 

0.522 

0.899 

0.896 

0.828 

0.930 

0.560 
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Two cases emerged from DOE 1 as the go-forward designs: jmp08 and optc. A comparison of 
some mechanical and geometry parameters is shown in Table 4.2. 3-4. Also shown in the table is 
another optimization design, opta, because it revealed an important aspect of integrating all the 
design disciplines. Although it had some highly desirable features in terms of stress and sweep, 
rotor-stator interaction noise results from the NASA - Acoustic code V072 showed the case to be 
unacceptable. (A comparison of the opta, optc and jmp08 results is shown with the DOE 2 
results, Section 4. 2.4.4). Results from DOE 1 go-forward cases; jmp08, opta, and optc showed 
that at the conclusion of DOE 1, the aerodynamic performance was below the design goal, and 
would have to be increased in the other DOE’s planned. Although the results of these cases 
looked promising, additional blade effective sweep was required. 

Table 4.2.3-4. DOE 1 Go-Forward Case Comparison 



Desired 

OPTA 

OPTC 

JMP08 

Optimization 


Acoustic 

Performance 

Aero 

Performance 

DOE #1 

Stresses 

<85 Ksi 

93 Ksi 

73 Ksi 

120 Ksi 

Defln. 

<0.4 in. 

0.35 in. 

0.47 in. 

0.55 in. 

YCG* Tip 


1.475 in. 

1.774 in. 

1.343 in. 

YCG* Mid 


-0.443 in. 

-0.152 in. 

-0.446 in. 

XCG** Tip 

Fwd (-) 

-1.5 in. 

-1.871 in. 

-1.871 in. 

XCG** Mid 

Fwd (-) 

40.281 in. 

-0.333 in. 

-0.279 in. 

OBA075 

> l-34°l 

i 

CD 

80 

o 

-25° 

-22° 

OBAIOO 

> 1-41 °l 

Ui 

o 

o 

-49° 

-41° 


* Positive Toward Pressure Side 

** Positive Toward Aft 


Late in the design process (after completion of the rotor design and half of the stator design 
effort), it was discovered that the V072 input parameter ACLS, the relative flow angle at the 
stator leading edge, was being input incorrectly. Time allowed for recalculation of the baseline 
engine calibration, the original rotor DOE 1, and all pertinent stator DOE cases. Although 
incorrect predictions were used to direct the rotor design process, the general conclusions derived 
from these analyses were supported by the corrected predictions. 
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V072 predictions were generated for the rotor DOE 1 cases 1-9 and the optimized cases jmp08 
and optc. All the rotor DOE cases used the baseline stator design. The inlet and aft duct power 
levels for 2 and 3 times the bladepass frequency are shown in Figure 4.2. 3-9 as deltas from the 
baseline predicted power levels. 



1 23456789 jmp8 optc 

Case 

Figure 4.2.3-9. V072 Predictions for Cases in Rotor DOE 1, All Use Baseline Stator 
(1*BPF cutoff) 
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4.2.4 Design of Experiment 2 

4.2.4. 1 Purpose 

The second DOE was created to minimize flutter. Although the goal was to use aeroelastic 
programs (FREPS, UNSFLO, and TURBO- AE) to calculate the aerodynamic damping at various 
points on the performance map, unfortunately the codes either had problems running or were 
considered not well enough calibrated to use in the DOE. Empirical data was used as the 
aeroelastic quality factors and direct the blade design. 

4.2.4.2 Description 

The second DOE consisted of 7 factors (position of max thickness at hub, mid, and tip; Tmax/C 
at hub, mid, and tip; and Ycg-tip) with 2 levels each. The 8 cases are described in Table 4.2.4-1. 
Although DOE 1 had shown that Ycg had little impact on aeroelastic empirical data, it was 
included to alleviate potential weight increase produced by thickness increases required to reduce 
stress and not flutter. Values of Ycg tip were selected based on DOE 1 results. Because the 
DOE 1 go-forward case optc was designed slightly after jmp08, and the need to initiate DOE 2; 
jmp08 was used as the QHSF baseline in DOE 2 even though the peak stresses and maximum tip 
deflection were above desired levels. Ranges for the other values in the table were selected 
based on previous analyses using the ANSYS optimization process completed to determine 
whether the blade design should be a 1/rev or a 2/rev, (Section 4.2.1). Spanwise thickness and 
position of max thickness for these analyses, compared to the baseline, are shown in Figures 
4.2.4- 1 and 4.2.4-2. 
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Chord Fraction (0=LE) 


Table 4.2.4-1. DOE 2 Aeroelastic Analysis X-Factors 


Case No. 

Tmax Loc 

Tmax Loc 

Tmax Loc 

Tmax/c 

Tmax/c 

Tmax/c 

Y eg* 

Hub 

Mid 

Tip 

Hub 

Mid 

Tip 

Tip 


1 

0.25 

0.30 

0.50 

0.04 

0.065 

0.030 

1.80 

2 

0.25 

0.30 

0.70 

0.08 

0.030 

0.015 

1.80 

3 

0.25 

0.65 

0.50 

0.08 

0.030 

0.030 

1.00 

4 

0.25 

0.65 

0.70 

0.04 

0.065 

0.015 

1.00 

5 

0.60 

0.30 

0.50 

0.08 

0.065 

0.015 

1.00 

6 

0.60 

0.30 

0.70 

0.04 

0.030 

0.030 

1.00 

7 

0.60 

0.65 

0.50 

0.04 

0.030 

0.015 

1.80 

8 

0.60 

0.65 

0.70 

0.08 

0.065 

0.030 

1.80 


*Y-cg are in inches and relative to Baseline engine size (30.7 in Dia.) positive in 
direction of rotation 


• DOE Levels 

0.8 



Figure 4.2.4-1. Levels Selected for Position of Maximum Thickness 
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• DOE Levels 



Figure 4.2.4-2. Levels Selected for Maximum Thickness 


4.2.4.3 Analysis 

The DOE 2 cases were analyzed using the same techniques and same response variables (Y- 
factors) as DOE 1, shown in Table 4.2. 3-2. 


4.2.4.4 Results 

For consistency, the same response variables are shown as presented in Section 4.2.4.4. Included 
also in the figures are the go-forward cases from DOE 1; jmp08, opta, and optc. Again, it is the 
sensitivities of the design and response variables which were the objective of the DOE. Results 
from DOE 2, as well as jmp08, opta, and optc are shown in Figures 4.2.4-3 through 4.2.4-6. 
Acoustic response variables for interaction noise and MPT’s are shown in Figure 4.2.4-3. Rotor- 
stator interaction noise was again calculated using V072 at 55.9% Nlc, and was compared 
relative to the baseline such that a negative value was desired. The figure shows why opta, 
which looked good based on the other discipline design criteria, was not used as a go-forward 
design. The reason why it had such high interaction noise levels was that in an effort to create 
more effective sweep the blade was allowed to sweep aft at midspan prior to sweeping forward. 
The reduction in rotor-stator axial spacing caused the increased noise. The figure shows the large 
amount of effective sweep being considered relative to the aft-swept baseline rotor. 
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1 2 3 4 5 6 7 8 JMP08 OPT A OPT C 

Case No. 



Case No. 

Figure 4.2.4-3. Result of the Acoustic Assessment of the Configurations for 
Rotor DOE 2. 




Case No. 

1 2 3 4 5 6 7 8 JMP08 OPT A OPT C 



Figure 4.2.4-4. Results of the Aerodynamic Assessment of the Configurations 

for Rotor DOE 2. 
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Case No. 



Figure 4.2A-5. Results of the Aeroelastic Assessment of the Configurations for 

Rotor DOE 2. 




Figure 4.2.4-6. Results of the Mechanical Assessment of the Configurations for 
Rotor DOE 2. 
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A new methodology was required to complete rotor DOE 2. Because the blade geometry was 
created by stacking airfoil section centers of gravity (eg) on a radial line from the engine center 
line through the hub section eg, any change to the hub section maximum thickness (tmax), or 
position of tmax would cause the hub eg to shift and the airfoil geometry to change - in 
particular, the leading edge (LE) metal angle. Because the LE metal angle was changing between 
the cases, the aerodynamic quality characteristics (flow, efficiency, etc.) could not be directly 
related to the x-factors. A long and unsuccessful attempt was made to stack the optc airfoil from 
DOE 1 on the LE which would, if successful, allow direct comparisons between the x-factors and 
y-factors of DOE 2. Although this technique was successful for the stator, the large hub slope of 
the rotor proved too much for AXCAPS and back-to-back geometry generation cases could not 
be duplicated. 

In order to solve this problem, a technique was developed which generated the blade geometry in 
an ANSYS script using the airfoil generation routines from AXCAPS. This procedure forced the 
airfoil section mean camber lines to be the same as optc. Ranges were determined for the x- 
factors, and the internal ANSYS optimization procedure was used to minimize flutter. This was 
essentially the same optimization procedure described in the process flowchart (Figure 4.2. 3-4) 
and used in DOE 1 and the first part of DOE 2, except AXCAPS was not used to generate the 
airfoil. Go-forward cases, based on the mechanical and empirical aeroelastic criteria, were then 
analyzed using DAWES to obtain the aerodynamic quality characteristics. The final go-forward 
case from DOE 2 was designated optm. Results of optm relative to the two go-forward cases 
from DOE 1 (jmp08 and optc) are shown in Figure 4.2.4-7. The figures show that relative to the 
baseline, substantial improvements in empirical flutter parameters were achieved while reducing 
weight; without any additional flow reduction. 
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Figure 4.2.4-7. 


Rotor summary from DOE 2 


4.2.5 Design of Experiment 3 
4.2.5. 1 Purpose 

The purpose of DOE 3 was to increase flow and achieve desirable efficiency characteristics by 
optimizing the incidence. 


4.2.5.2 Description 

The base level used for DOE 3 was optm. In order to optimize the incidence across the radial 
span, five spans were chosen (0%, 25%, 50%, 75%, and 95%) with three levels at each span as 
shown in Table 4.2.5- 1. The tip (100% span) was not included because of an inlet boundary 
layer which was being modeled in the analysis; however, the tip was adjusted to remain 
consistent with the 95% span level. The magnitude in incidence was also reduced from hub to 
95% span to satisfy an additional flutter criteria involving excessive positive tip incidence. The 
minimum fractional factorial orthogonal matrix for 5 factors with 3 levels produced 27 different 
geometry’s; and, because four or more runs were required for each geometry to determine the 
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speedline shapes, it resulted in over 108 runs. Because of the time required to run such a large 
number of cases, DENTON (3-D inviscid analysis) was again reviewed as a possible analysis 
tool. 


Table 4.2.5-1. DOE 3 Design Point Rotor Performance Analysis 


Incidence Values in Degrees 
(Values represent differences from base case - OPTM) 


Case # 

0 % Span 

25 % 

50 % 

75 % 

95 % 

1 

- 4.0 

- 3.5 

- 3.0 

- 2.5 

- 2.0 

2 

- 4.0 

- 3.5 

0.0 

2.5 

2.0 

3 

- 4.0 

- 3.5 

3.0 

0.0 

0.0 

4 

- 4.0 

0.0 

- 3.0 

2.5 

2.0 

5 

- 4.0 

0.0 

0.0 

0.0 

0.0 

6 

- 4.0 

0.0 

3.0 

- 2.5 

- 2.0 

7 

- 4.0 

3.5 

- 3.0 

0.0 

0.0 

8 

- 4.0 

3.5 

0.0 

- 2.5 

- 2.0 

9 

- 4.0 

3.5 

3.0 

2.5 

2.0 

10 

0.0 

- 3.5 

- 3.0 

2.5 

0.0 

11 

0.0 

- 3.5 

0.0 

0.0 

- 2.0 

12 

0.0 

- 3.5 

3.0 

- 2.5 

2.0 

13 

0.0 

0.0 

- 3.0 

0.0 

- 2.0 

14 

0.0 

0.0 

0.0 

- 2.5 

2.0 

15 

0.0 

0.0 

3.0 

2.5 

0.0 

16 

0.0 

3.5 

- 3.0 

- 2.5 

2.0 

17 

0.0 

3.5 

0.0 

2.5 

0.0 

18 

0.0 

3.5 

3.0 

0.0 

- 2.0 

19 

4.0 

- 3.5 

- 3.0 

0.0 

2.0 

20 

4.0 

- 3.5 

0.0 

- 2.5 

0.0 

21 

4.0 

- 3.5 

3.0 

2.5 

- 2.0 

22 

4.0 

0.0 

- 3.0 

- 2.5 

0.0 

23 

4.0 

0.0 

0.0 

2.5 

- 2.0 

24 

4.0 

0.0 

3.0 

0.0 

2.0 

25 

4.0 

3.5 

- 3.0 

2.5 

- 2.0 

26 

4.0 

3.5 

0.0 

0.0 

2.0 

27 

4.0 

3.5 

3.0 

- 2.5 

0.0 


4.2.5.3 Analysis 

DENTON was run with the optc geometry and speedlines were compared with those previously 
completed using DAWES. Although the total pressure and efficiency levels were considerably 
different (as expected), choke flow and speedline shapes of efficiency and pressure compared 
reasonably well. Therefore DENTON was used to complete DOE 3. After completing the 
analysis, the JMP program and EXCEL Solver were used to predict go-forward cases which 
would achieve the desirable results. The y-factors were choke flow, a flow delta between deep 
choke and design-point, and a qualitative assessment of the efficiency speed line shape. 
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4.2.5.4 Results 


Statistical correlations of the data were not very good, but several go-forward cases were 
completed. The final go-forward case from DOE 3 was designated optm3 and is shown 
compared to three representative cases from the DOE: Cases 1, 9, and 19 in Figure 4.2.5- 1. The 
figure shows that, in general, cases with large amounts of choke flow (Case 9) had poor 
efficiency shapes that reduced in level from choke to stall, while those with lower amounts of 
choke flow (Case 1) had efficiency shapes that increased in level from choke to stall. The go- 
forward case optm3 was selected to provide maximum flow while still exhibiting good efficiency 
speedline shape characteristics. 




Wc (Ibm/s) (Baseline engine flow size) Wc (Ibm/s) (Baseline engine flow size) 


Figure 4.2.5-1. Representative Results from DOE 3 Compared to the Go- 
Forward Case OPTM3. 
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4.2.6 Design of Experiment 4 


4.2.6. 1 Purpose 

The purpose of DOE 4 was to increase design point flow and efficiency by optimizing the blade 
metal angle (beta) distributions along the mean camber. 


4.2.6.2 Description 

A fractional factorial L-16 was used incorporating fifteen x-factors with two levels each. 
(Although 3 levels for each x-factor would have been preferred, completion time associated with 
using DAWES prevented it.) Go-forward case optm3 (from DOE 3) was used as the base level. 
DAWES analyses performed on optm3 showed that efficiency levels could be increased across 
the span and therefore dictated the complete spanwise range for DOE 4 design factors. A 
description of the x and y-factors selected are included in Table 4. 2.6-1 along with the matrix 
description in Table 4.2.6-2. The x-factors are shown graphically in Figure 4.2.6-1 as ranges 
from the base configuration (optm3). For the blade metal angles on streamlines between those 
being modified, the betas were linearly interpolated for a smooth transition. 
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Table 4.2.6- 1. DOE 4 Design Point Rotor Performance Analysis 
Quality Characteristics 


(Values represent differences from base case - OPTM3) 


X-factors 

+/- 3.0 deg. 
+/- 3.0 deg. 
+/- 3.0 deg. 
+/- 2.5 deg. 
+/- 2.5 deg. 
+/- 2.5 deg. 
+/- 2.5 deg. 
+/- 2.0 deg. 
0, -20.5 deg. 
+/- 2.0 deg. 
+/- 2.0 deg. 
15%, 25% 
40%, 60% 
10%, 25% 
40%, 60% 


Y-f actors 

Design point flow 
Peak efficiency 

Design point efficiency radial profile 
(10, 30, 50, 70, 90% radial spans) 
Aerodynamic loadings 

(10, 30, 50, 70, 90% radial spans) 
Suction surface passage shock location 
(50, 70, 90% radial spans) 

Surge margin assessment 
Part-speed efficiency 


1 Beta SL1 LE 

2 Beta SL1 25% m 

3 Beta SL1 50% m 

4 Beta SL12 LE 

5 Beta SL1 2 25% m 

6 Beta SL1 2 50% m 

7 Beta SL1 2 75% m 

8 Beta SL24 LE 

9 Slope of Beta to x-f 1 4 loc. 

1 0 Beta SL24 50% m 

1 1 Beta SL24 75% m 

1 2 Beta location of x-factor 5 

1 3 Beta location of x-factor 6 

1 4 Beta location of x-factor 9 

1 5 Beta location of x-factor 1 0 


Table 4.2.6-2. DOE 4 Design Point Rotor Performance Analysis 

X-Factor Levels 


(Values represent differences from base case - OPTM3) 


Case # 

XI 

X2 

X3 

X4 

X5 

X6 

X7 

X8 

X9 

XI 0 

XII 

XI 2 

XI 3 

XI 4 

XI 5 

1 

- 3.0 

- 3.0 

- 3.0 

- 2.5 

2.5 

- 2.5 

- 2.5 

2.0 

- 20.5064 

2.0 

2.0 

4.0 

13.0 

6.0 

13.0 

2 

- 3.0 

- 3.0 

- 3.0 

2.5 

- 2.5 

2.5 

2.5 

- 2.0 

0.0 

- 2.0 

- 2.0 

4.0 

13.0 

6.0 

13.0 

3 

- 3.0 

- 3.0 

3.0 


- 2.5 

2.5 

2.5 

- 2.0 

- 20.5064 

2.0 

2.0 

6.0 

9.0 

3.0 

13.0 

4 

- 3.0 

- 3.0 

3.0 


2.5 

- 2.5 

- 2.5 

2.0 

0.0 

- 2.0 

- 2.0 

6.0 

9.0 

3.0 

13.0 

5 

- 3.0 

3.0 

- 3.0 

- 2.5 

- 2.5 

2.5 

- 2.5 

2.0 

0.0 

- 2.0 

2.0 

6.0 

9.0 

6.0 

9.0 

6 

- 3.0 

3.0 

- 3.0 

2.5 

2.5 

- 2.5 

2.5 

- 2.0 

- 20.5064 

2.0 

- 2.0 

6.0 

9.0 

6.0 

9.0 

7 

- 3.0 

3.0 

3.0 

- 2.5 

2.5 

- 2.5 

2.5 

- 2.0 

0.0 

- 2.0 

2.0 

4.0 

13.0 

3.0 

9.0 

8 

- 3.0 

3.0 

3.0 

2.5 

- 2.5 

2.5 

- 2.5 

2.0 

- 20.5064 

2.0 

- 2.0 

4.0 

13.0 

3.0 

9.0 

9 

3.0 

- 3.0 

- 3.0 

- 2.5 

- 2.5 

- 2.5 

2.5 

2.0 

0.0 

2.0 

- 2.0 

6.0 

13.0 

3.0 

9.0 

10 

3.0 

- 3.0 

- 3.0 

2.5 

2.5 

2.5 

- 2.5 

- 2.0 

- 20.5064 

- 2.0 

2.0 

6.0 

13.0 

3.0 

9.0 

11 

3.0 

- 3.0 

3.0 

- 2.5 

2.5 

2.5 

- 2.5 

- 2.0 

0.0 

2.0 

- 2.0 

4.0 

9.0 

6.0 

9.0 

12 

3.0 

- 3.0 

3.0 

2.5 

- 2.5 

- 2.5 

2.5 

2.0 

- 20.5064 

- 2.0 

2.0 

4.0 

9.0 

6.0 

9.0 

13 

3.0 

3.0 

- 3.0 

- 2.5 

2.5 

2.5 

2.5 

2.0 

- 20.5064 

- 2.0 

- 2.0 

4.0 

9.0 

3.0 

13.0 

14 

3.0 

3.0 

- 3.0 

2.5 

- 2.5 

- 2.5 

- 2.5 

- 2.0 

0.0 

2.0 

2.0 

4.0 

9.0 

3.0 

13.0 

15 

3.0 

3.0 

3.0 

- 2.5 

- 2.5 

- 2.5 

- 2.5 

- 2.0 

- 20.5064 

- 2.0 

- 2.0 

6.0 

13.0 

6.0 

13.0 

16 

3.0 

3.0 

3.0 

2.5 

2.5 

2.5 

2.5 

2.0 

0.0 

2.0 

2.0 

6.0 

13.0 

6.0 

13.0 
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Figure 4.2.6-1. Graphical Description of the X-Factors Used in Rotor DOE 4. 


4.2.6.3 Analysis 

DAWES results were analyzed and compared for the 64 runs (16 cases with 4 different 
backpressures) at the design speed. Y-factors of flow, efficiency, pressure ratio, and work were 
calculated and plotted. Span wise efficiency profiles were compared and tabulated at 10, 30, 50, 
70, and 90 percent span. Aerodynamic blade loadings were also compared at the same spans and 
were qualitatively assigned values for the DOE. The loadings were judged on suction surface 


NAS A/CR— 2003-2 1 2369 


37 



peak Mach number, suction surface diffusion rate, leading edge assessment (choke/stall), and 
pressure surface Mach number distribution. Passage shock locations (on the suction surface) 
were calculated and tabulated at 50, 70, and 90 percent span. Once the y-factors had been 
tabulated, JMP and EXCEL were used to create go-forward cases. 

One potential problem in using DOE’s with two levels is that it assumes that the solution is 
linear, which is typically not the case. This results in extrapolated values with excessive ranges 
which would not be feasible. To address this, constraints were imposed on the amount that the 
design factors could permeate. This forced the optimization procedure to find a feasible solution 
(all design criteria met) within the specified range. 


4.2.6.4 Results 

The tabulated results are shown in Table 4.2.6-3. A representative sample of the cases (low flow, 
high flow, and max efficiency) are shown at the design point relative to the baseline rotor 
DAWES analysis in Figure 4.2. 6-2. 


Table 4.2.6-3. DOE 4 Tabulated Results 


Case 

Wc 17.0 

Eff pk 

Eff 10% 

Eff 30% 

Eff 50% 

Eff 70% 

Eff 90% 

L 10% 

L30% 

L 50% 

L 70% 

L 90% 

S L 50% 

S L 70% 

S L 90% 

SM 

1 

183.78 

0.8914 

0.9797 

0.9415 

0.8992 

0.8784 

0.8136 

13 

9 

2 

8 

2 

0.247 

0.343 

0.649 

-0.0215 

2 

188.24 

0.8931 

0.9767 

0.9279 

0.8829 

0.8831 

0.8451 

12 

14 

11 

10 

13 

0.246 

0.170 

0.701 

-0.0094 

3 

185.40 

0.8772 

0.9762 

0.9055 

0.8425 

0.8595 

0.8783 

16 

16 

16 

16 

8 

0.166 

0.189 

0.595 

0.0000 

4 

180.95 

0.8832 

0.9819 

0.9474 

0.8956 

0.8556 

0.6982 

14 

6 

6 

13 

10 

0.291 

0.418 

0.596 

-0.0074 

5 

183.33 

0.8845 

0.9790 

0.9225 

0.8844 

0.8734 

0.8294 

15 

15 

13 

6 

12 

0.185 

0.343 

0.650 

-0.0215 

6 

188.68 

0.9040 

0.9802 

0.9493 

0.8948 

0.8916 

0.8405 

11 

3 

14 

2 

16 

0.339 

0.210 

0.700 

-0.0128 

7 

190.02 

0.8994 

0.9803 

0.9458 

0.8983 

0.8931 

0.8603 

10 

5 

3 

1 

5 

0.269 

0.296 

0.699 

0.0000 

8 

186.04 

0.8965 

0.9795 

0.9306 

0.8830 

0.8855 

0.8616 

9 

11 

12 

5 

7 

0.225 

0.189 

0.675 

-0.0322 

9 

184.90 

0.8859 

0.9775 

0.9300 

0.8840 

0.8775 

0.8416 

4 

12 

4 

4 

3 

0.205 

0.295 

0.649 

-0.0196 

10 

188.76 

0.9072 

0.9790 

0.9462 

0.8925 

0.8836 

0.8589 

1 

7 

15 

12 

9 

0.316 

0.210 

0.700 

-0.0120 

11 

186.55 

0.8984 

0.9801 

0.9456 

0.8944 

0.8689 

0.8431 

5 

2 

10 

15 

15 

0.248 

0.134 

0.700 

-0.0149 

12 

187.91 

0.8909 

0.9773 

0.9355 

0.8830 

0.8771 

0.8479 

3 

13 

7 

14 

1 

0.291 

0.319 

0.699 

-0.0139 

13 

183.81 

0.8995 

0.9786 

0.9518 

0.9034 

0.8795 

0.8649 

8 

1 

1 

11 

4 

0.248 

0.418 

0.675 

-0.0314 

14 

189.80 

0.8982 

0.9803 

0.9419 

0.8897 

0.8926 

0.8353 

7 

8 

8 

3 

14 

0.292 

0.189 

0.722 

-0.0098 

15 

194.05 

0.8970 

0.9791 

0.9380 

0.8848 

0.8925 

0.8766 

6 

10 

5 

7 

6 

0.247 

0.250 

0.674 

0.0000 

16 

177.01 

0.8938 

0.9797 

0.9458 

0.8957 

0.8612 

0.6605 

2 

4 

9 

9 

11 

0.269 

0.446 

0.700 

-0.0087 
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Figure 4.2.6-2. Representative Results from Rotor DOE 4 Compared to 
Baseline Rotor 

As Table 4.2.6-3 shows, the calculated efficiency and the qualitative value given to the 
aerodynamic loading, for the same spanwise location, varied considerably in some instances. 
Therefore, go-forward cases were created based on the other y-factors combined with either the 
efficiency at radial spans or the loadings. DAWES results for the go-forward cases at the design 
point are shown in Figure 4. 2.6-3 compared to the starting case from DOE 3 (optm3) and the 
performance objective (baseline rotor). The initial two cases (gfOl, and gf02) were completed 
using EXCEL and were predicted to match the design flow objective; however, DAWES results 
showed insufficient flow. Y-factors from these cases were added to the model and two new go- 
forward cases were created, gf03 and gf04. The four cases were also analyzed at 85% Nlc to 
verify good part-speed performance with results shown in Figure 4. 2.6-4. The results indicated a 
reduction in part-speed efficiency with increased choke flow at 100% Nlc. In order to best 
match the baseline rotor performance goals, and provide the best back-to-back test for acoustic 
comparisons, go-forward case 3 (gf03) was selected as the final design (details reported in 
Section 5.0). 
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Rotor PR Rotor PR 



182 184 186 188 190 192 194 196 

Wc (Ibm/s) (Baseline engine flow size) 



Wc (Ibm/s) (Baseline engine flow size) 


Figure 4.2.6-3. Effect of Beta Distribution on Flow and Efficiency on DOE 4 Go- 
Forward Cases at 100% Nlc. 




Figure 4.2.6-4. Effect of Beta Distribution on Flow and Efficiency on DOE 4 Go- 
Forward Cases at 85% Nlc. 
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4.2.7 Bird Ingestion Analysis 


A preliminary Foreign-Object-Damage (FOD) assessment performed on the QHSF blade 
indicated leading-edge thickness parameters to be well within successful AE experience. This 
analysis is comparative in nature and considers only spanwise geometric characteristics, LE 
thickness distribution, blade count, metal angles, bird weight, as well as other parameters. The 
computed damage tolerance factors are compared with a database consisting of very successful, 
marginal, and poor designs. This preliminary analysis however, can not account for the high 
degree of forward sweep and tangential lean built into the QHSF, which will be accounted for 
using NOSAPM (Section 5. 1.4.6). 
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4.3 Stator 


4.3.1 Design Ground Rules and Goals 

The stator design emphasizes the reduction of rotor wake interaction noise while meeting 
aerodynamic, mechanical, and producibility goals. The ground rules required: 

1. ) Use the baseline vane count (52 vanes). The acoustic effects of varying vane count would 

confuse the analysis of the effects of vane design. 

2. ) Don’t sweep the stator aft or move it downstream relative to the baseline stator. The fan will 

be able to replace the baseline fan without a change in engine dimensions. 

3. ) Use the baseline front frame unchanged. 

4. ) Use the baseline stator solidity. The inflow from the AE/QHSF was similar to the baseline 

fan so that the baseline solidity would provide good loadings. 

5. ) Control vane stacking by specifying the shape of the leading edge. Since this is what affects 

wake trace speed, it seemed simpler to specify the shape of the leading edge. 

6. ) Design for composite construction. This requirement puts some restrictions on leading and 

trailing edge radius and maximum thickness. 


The stator was to reduce rotor wake interaction noise to the greatest extent possible while still 
having acceptable aerodynamic loadings, losses similar to the baseline stator, and acceptable 
stress, frequency margin, flutter margin, and composite producibility. 


4.3.2 Strategy 

The vane stacking would be set first since this was important acoustically and since it was likely 
to be radical enough that it would be good to find out up front about any aerodynamic or 
mechanical problems it caused. DOE 1 would be used to experiment with different vane stacks. 
As originally planned, this DOE varied vane stacking tangentially and axially. Axial stacking 
changes were outside the ground rules, but knowledge of these trends was desirable anyway. 
Acoustic, aerodynamic, and mechanical quality factors would be collected and their trends with 
vane stacking would be used to design the final vane stack. The vane incidence, vane section 
mean camber line angle chordwise distributions, deviation, chord, and thicknesses for this DOE 
would come from the baseline vane in the beliefs that it was a good starting point. Aerodynamic 
quality characteristics would be calculated by the DAWES 3D viscous CFD code. 
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Once the stacking had been decided, the second vane DOE would be used to optimize vane 
incidence and mean camber line angle chordwise distributions for good aerodynamic 
performance with this stacking. 

The vane trailing edge angle would then be set to produce axial (zero swirl) stator exit flow. 
This would not require a DOE. The vane trailing angle would be set to remove any residual 
DAWES -predicted exit swirl. This angle usually produces exit swirl acceptably close to axial in 
just a few iterations. Loadings and loss are then re-checked to make sure they’re still acceptable 
in spite of any increased flow turning required to produce axial flow. 

The vane leading and trailing edge radius and maximum thicknesses would then be adjusted to 
meet composite vane producibility criteria. Composite vanes require control of the vane 
thickness as measured normal to the vane suction and pressure surfaces, as opposed to the design 
tools which specify thickness in the tangential direction. The tangential thickness varies from the 
normal thickness by the cosine of the lean angle. Since the lean angle everywhere but the leading 
edge is affected by both the vane stacking and the vane section stagger angles, the thickness 
modification is performed after the stagger angles have been determined. After thickness has 
been adjusted, the loadings, loss, and exit swirl are re-checked. 

The composite construction provides flexibility to produce a vane with combinations of high 
axial and tangential leans with minimum weight penalty. The current production vanes on the 
baseline stator have symmetric lay-up or taped construction. Recently, a braided vane has been 
developed at AE as a part of the ENCORE (Engine Cost Reduction) program. In addition, the 
braided vane has better FOD characteristics than the taped vane. However, it has lower flutter 
margins due to decreased torsional rigidity. The QHSF exit vane will be a braided composite 
vane. A preliminary layout of the vane is shown in Figure 4. 3. 2-1. 

The AE in-house design tool CADS (Composite Airfoil Design System) will be used for the 
composite vane design. As shown in the flowchart in Figure 4. 3. 2-2, CADS uses the airfoil 
geometry from the aero bank file to generate an ANSYS input (prep7) file. This system was used 
to design the baseline stator. 
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Band 




Figure 4.3.2-2. The AE Composite Stator Design System will Design a 
Composite Vane that Will Incorporate the Optimum 
Aerodynamic and Acoustic Features. 
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4.3.3 Stacking DOE 


The vanes used in the stacking DOE were designed to operate behind the “optc” fan rotor. 
Initially, an LI 8 Taguchi matrix was chosen with 2 levels for axial sweep at the vane midspan 
and tip and 3 levels for the tangential stacking at the midspan and tip. The 2 sweep levels 
represented the baseline vane axial sweep and a 14 degree increase. This DOE couldn’t be 
completed because the great variation in rotor-stator axial spacing between the hub and the 
shroud created convergence problems in the axisymmetric flow solver (AXCAPS). The matrix 
was modified so the sweep now varied between the baseline sweep and a radial leading edge, 
which would still provide information about the effects of sweep. Work on the modified DOE 
was abandoned due to schedule constraints and since decreasing axial sweep had noticeable 
acoustic and aerodynamic penalties. 

The baseline axial sweep was then chosen for the vane and only the effect of tangential lean and 
bow was studied. Tangential stacking was controlled at the midspan and the shroud with three 
levels at each location representing lean with the direction of rotor rotation, no lean, and lean 
against the direction of rotor rotation. Since this only produced 9 different vanes, a full factorial 
was used to examine all 9 vanes. 

In full factorial stator DOE 1 (FF 1), tip tangential displacement equivalent to 30 deg of lean 
from the radial were tried. This had acoustic benefits but vanes leaned 30 deg with the rotor 
rotation (i.e., with the suction surface leaned down facing the hub) showed serious flow 
separation near the hub in DAWES. Since the 30 deg lean was infeasible from a performance 
standpoint, a second full factorial, FF 2, was used which had vane section tangential tip 
displacements equivalent to 15 deg of lean from the radial. The leading edges are shown in Fig. 
4.3. 3-1 and tabulated in Table 4. 3. 3-1. This figure also shows rev_3, which was the stacking 
eventually chosen after vane iterations using the acoustic, aerodynamic, and mechanical results 
of FF 2. 
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Figure 4.3.3-1. Stator Stacking Full Factorial 2 Vane Leading Edges 


Table 4.3.3- 1. Stator Stacking Full Factorial 2 

CG Shifts in Baseline Engine Scale 


Case 

1 

Midspan Ycg (in) 
-0.879 

Shroud Ycg (in) 
-1.758 

15 deg lean with rotation 

2 

0.000 

0.000 

No lean 

3 

0.879 

1.758 

15 deg lean against rotation 

4 

-0.879 

0.000 

Bow with rotation 

5 

-0.879 

1.758 

Compound bow w/ and against rotation 

6 

0.000 

-1.758 

Compound lean with rotation 

7 

0.000 

1.758 

Compound lean against rotation 

8 

0.879 

-1.758 

Compound bow against and w/ rotation 

9 

0.879 

0.000 

Bow against rotation 


Ycg: +ve=shifted against direction of rotation 
0=unshifted 

-ve=shifted with direction of rotation 
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The quality characteristics used for the vane stacking iterations were: 

Acoustic: V072-predicted rotor wake-stator interaction sound power levels (SPL) at 2 and 3 
times blade passing frequency (BPF), radiated towards the fan inlet and exit, at a representative 
approach point at 55.9% corrected fan design speed were used. The stator inlet velocity triangles 
used in V072 were predicted by an off-design model of the optc rotor in AXCAPS. The acoustic 
results are shown in Figures 4.3. 3-2 through 4.3. 3-5 for inlet 2*BPF, inlet 3*BPF, exit 2*BPF, 
and exit 3*BPF respectively. These figures also show the acoustic predictions for rev_3, the 
final vane stack. 



Figure 4.3.3-2. Stator Stacking Full Factorial 2 Inlet 2*BPF SPL Predictions 
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Figure 4A.3-5. Stator Stacking Full Factorial 2 Exit 3*BPF SPL Predictions 


The acoustic results were difficult to optimize. Generally, any stacking that reduced one sound 
power level increased others. A reduction in inlet sound power level was favored over a 
reduction in exit because the exit radiated sound can be more easily attenuated with acoustic 
treatment. The inlet 2*BPF results favored leaning the midspan against rotation and the tip with 
rotation as in rev_3. 

Aerodynamic: Design point stator loss and stator exit swirl as estimated by DAWES were used 
as quality factors. Stator loadings were used as a qualitative guide. Figure 4.3. 3-6 shows the 
DAWES predicted loss for the 9 stators in FF 2 as well as for the rev_3 stator. 
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Figure 4.3.3-6. Stator Stacking Full Factorial 2 Stator Loss 


The loss profiles show that leaning the suction surface towards the endwalls should be avoided 
because of flow separation from the suction surface and high hub loss (Cases 1, 4, and 5) and 
shroud loss (Cases 5 and 7). This result was an important aerodynamic driver in the stacking 
iterations. 

The loss profiles also show a high-loss region from about 60% to 80% span for Cases 6, 8, and 
(to a lesser extent) 9. At first, this appeared to be a problem with stacks that lean with rotation 
over the outer 50% span but examination of the flow vectors in the vane-to-vane surface in this 
region showed flow separation from the pressure surface due to large negative incidence. Figure 

4.3. 3- 7 shows incidence calculated using DAWES leading edge swirl angles. Figures 4.3. 3-8 and 

4.3. 3- 9 show the vane leading edge geometric angles and the vane leading edge swirl angles that 
resulted in this incidence. 
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Figure 4.3.3-7. Stator Stacking Full Factorial 2 Stator Incidence 



Figure 4.3.3-8. Stator Stacking Full Factorial 2 Vane LE Angle 
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Figure 4.3.3-9. Stator Stacking Full Factorial 2 Stator LE Swirl Angle 

The vane leading edge geometric angles in Figure 4.3. 3-8 were calculated by AXCAPS. The 
intended stator incidence for FF 2 was the baseline design point stator incidence. AXCAPS was 
used to calculate the required vane leading edge angles to achieve this incidence given the rotor 
exit swirl angles from the OPTC rotor DAWES model. The varying vane stacks produced 
varying stator leading edge swirl profiles in AXCAPS since AXCAPS applies vane body radial 
forces to shift the flow radially, which changes the stator leading edge meridional velocity profile 
according to how the vane is stacked. Since all 9 cases have the same incidence, the variation in 
vane leading edge geometric angle in Figure 4.3.3-8 represents the variation in leading edge swirl 
angle in AXCAPS, which can be slightly over 5 degrees between 60% and 80% span. 

The air angles in Figure 4.3. 3-9 were obtained from the stator DAWES models. Each DAWES 
model used the exit swirl profile predicted by the OPTC rotor DAWES model at the design 
point. The swirl profile variation from case to case at the stator leading edge represents 
DAWES’ estimate of how vane stacking affects the radial distribution of flow upstream of the 
stator. In DAWES, the variation due to vane stacking is only 2-3 degrees. 

The high loss in this region results from setting the vane leading edge angles too high for these 
cases. This happened because the streamline code which was used to generate vane geometry 
over-estimated the effect of vane stacking on the stator leading edge swirl angles. Stator loss in 
this region was, therefore, not considered in the vane stacking evaluation. Figure 4.3. 3-6 shows 
that the final vane stack, rev_3, also has high loss. This issue was resolved in the next DOE 
where incidence would be optimized. 
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Figure 4.3.3-10 shows stator exit swirl estimated by DAWES. The ultimate intent of the design 
was to set this as close to axial (0 deg) as practical, which was done in the baseline vane design. 
The vane trailing edge geometric angles for this DOE were taken from the baseline vane and 
were not adjusted to reduce the exit swirl. Exit swirl was still used as a quality characteristic 
since stators with swirl already close to axial would require less trailing edge angle adjustment. 



Figure 4.3.3-10. Stator Stacking Full Factorial 2 Stator TE Swirl Angle 


Since all the vanes in this study have the same trailing edge geometric vane angle, Figure 4.3.3- 
10 shows the expected result that vanes with higher loss have higher deviation. Although exit 
swirl was examined in the design process, it was not used directly in the DOE since it is another 
way of representing loss. 

Vane displacements, static strains on pressure and suction surfaces, vibration and flutter margins 
were used in the DOE as mechanical inputs. In addition, the static stresses and strains in the 
vane foot were monitored for abnormalities. The strain and displacement limits were defined 
based on AE experience. As seen from Figures 4.3.3-11 and 4.3.3-12, there were no surprises 
with both strain levels and the maximum displacements. Further, the maximum strains and 
displacements were within acceptable limits for all the go-forward cases that emerged out of the 
stator DOE. 
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Deflection in Inches Micro Strain 


1400 



Figure 4.3.3-11 Stator Stacking Full Factorial Maximum Static Strain 
Predictions 
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Figure 4.3.3-12 Stator Stacking Full Factorial Maximum Deflection Predictions 
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REV3 


Satisfying the AE empirical flutter criteria for vanes conflicted with the aerodynamic and 
acoustic design. From an aerodynamic/acoustics perspective, a high degree of tangential bow is 
desirable. However, the addition of bow causes a significant change in the mode shape of the 
fundamental vibration mode. With small amounts of bow this mode is predominantly a flexure 
mode, while a larger degree of bow changes it into a torsion mode (see Figure 4.3.3-13). The 
problem arises because the flutter criteria based on reduced frequency, fc/V, is dependent on 
mode shape, with the torsion mode being much more difficult to satisfy. Fundamental mode 
flutter parameters were found to be at least 25% lower than AE recommended values for most of 
the go-forward designs (see Figure 4.3.3-14). 

A number of approaches were investigated to increase the frequency of the bowed stator vane, 
such as thickness increases, changing airfoil chord, attachment boundary conditions, etc. Table 
4.3. 3-2 summarizes the flutter predictions for all the cases investigated. Initial attempts to 
increase flutter margin by thickening the vanes resulted in rapidly increasing stator weight with 
some moderate aerodynamic loss increases predicted by DAWES. For example, the total weight 
of the case_6 stator (52 vanes) increased from 5.68 to 7.07 lb. to drive the torsion mode flutter 
from 0.89 to 0.96. Increasing the chord further increased the weight with negligible change to 
the flutter margins. The weight issue was considered a problem and the low flutter margin was 
not necessarily a problem. There was no data for stators that implied flutter was guaranteed at 
this level of the flutter parameter, it was simply below previous experience. A compromise was 
reached where the vanes would be left at the rev_3 flutter parameter level and the vanes would be 
strain gauged in the rig to watch for flutter. Vane hub and shroud attachment modifications have 
been identified analytically that will increase flutter margin to baseline levels if vane flutter is 
seen in the rig. The modification is to provide two retention bands instead of the conventional 
one-band retention scheme used in the previous designs. The approach was tested on a number 
of cases and was found to consistently increase the flutter parameter by about 10 to 13% (see 
Table 4.3. 3-2). This approach helps keep stator weight under control and provides new 
knowledge about the limits of vane flutter parameter while providing an acceptable alternative if 
flutter actually occurs. 
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Flutter Parameter 
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(b) Rev_3 



Figure 4.3.3-13 Effect of Tangential Bow on the Mode Shape of Fundamental 
Vibration Mode 
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Figure 4.3.3-14. Stator Stacking Full Factorial Fundamental Mode Flutter 
Predictions 
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Table 4.3.3-2. Summary of Fundamental Modal Frequencies for Stator Vane 


DOE Case 


Mode 1 


Total Wt. (lb.) 

Comment 


Frequency 

(N/rev) 

Type* 

Flutter 

Parameter 



Case_l 

3.31 

B 

0.96 

5.64 


Case_2 

3.28 

B 

0.93 

5.60 


Case_3 

2.96 

B 

0.81 

5.63 


Case_4 

2.39 

B 

0.67 

5.63 


Case_5 

2.26 

B/T 

0.62 

5.56 


Case_6 

3.01 

T 

0.89 

5.68 


Case_6_Bnd 

3.29 

T 

0.97 

5.68 

Two Bands 

Case_6a 

3.27 

T 

0.96 

7.07 

Inc. Thickness 

Case_7 

2.37 

B 

0.66 

5.62 


Case_8 

2.14 

T 

0.64 

5.73 


Case_8c 

2.50 

T 

0.76 

8.45 

Inc. Thickness 

Case_9 

2.56 

T 

0.74 

5.65 


ESP_2 

2.49 

T 

0.74 

5.71 


ESP_2a 

2.83 

T 

0.85 

7.60 

Inc. Thickness 

ESP_2b 

2.82 

T 

0.84 

8.00 

Inc. Thk & Chord 

ESP_3 

2.18 

T 

0.65 

5.72 


REV_1 

2.62 

T 

0.77 

5.68 


REV_l_Bnd 

2.94 

T 

0.87 

5.68 

Two Bands 

REV_2 

2.72 

T 

0.80 

5.67 


REV_2_Bnd 

3.04 

T 

0.89 

5.67 

Two Bands 

REV_3 

2.75 

T 

0.81 

5.69 


REV_3_Bnd 

3.11 

T 

0.92 

5.69 

Two Bands 


* Type: B ^ Bending Mode; T ^ Torsion Mode 
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4.3.4 Incidence DOE 


The next step has to set the incidence and mean camber line angle distribution. The vane mean 
camber line angle distributions were taken from the baseline vane. Although the mean camber 
line angle distribution optimization would render reduced stator loss and improved stator 
loadings, the stator losses predicted by DAWES for the rev_3 stator were already close to those 
predicted by DAWES for the baseline stator. Since the final rotor design (gf03) was available, 
the AXCAPS model was modified to use this rotor with its DAWES -predicted exit velocity 
triangles for the remaining stator design work. 

The rev_3 stator behind the gf03 rotor was used as the base case for the incidence study. The 
vane leading edge angles were set to produce the baseline design point incidence with the stator 
leading edge swirl angles from the new gf03 rotor. A 9-case full factorial DOE was performed 
where vane leading edge angle was varied over three levels (-3 deg, 0 deg, and +3 deg) from the 
baseline at the vane hub, midspan, and tip. The vane leading edge angle was specified directly 
rather than allowing AXCAPS to try to calculate it from incidence. Figure 4.3. 4-1 shows the 
vane leading edge geometric angles for the 9 cases (Case 5 is the baseline case) and for the rev_3 
vane. Figure 4. 3.4-1 also shows the vane leading edge angles for inc_l, which is the leading 
edge chosen from the results of this study. Figure 4.3. 4-2 shows the resulting incidence profiles 
using DAWES stator leading edge swirl angles. 

The effects of varying the vane leading edge angle on the mechanical and acoustic quality factors 
were insignificant so the leading edge angle was chosen based on stator loss and loading. Exit 
swirl was examined, but for a constant vane trailing edge angle, the deviation largely followed 
the loss as in the stacking study. Figures 4.3.4-3 and 4. 3.4-4 show the stator loss and exit swirl 
estimated by DAWES. 
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Figure 4.3.4-1. Stator Incidence Full Factorial Incidence 



incidence (deg) from DAWES LE swirl 

Figure 4.3.4-2. Stator Incidence FF Stator Incidence 
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Figure 4.3.4-3. Stator Incidence FF Stator Loss 



Figure 4.3.4-4. Stator Incidence FF Stator TE Swirl Angle 
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The loss was reduced from or equal to rev_3 over much of the span for all 9 cases, indicating that 
the biggest benefit of the incidence full factorial was due to correcting the way in which 
incidence was set. The leading edge angle was set by looking at curves of loss vs. leading edge 
angle and at vane loadings at different spanwise locations. Generally, leading edge angles that 
improve loss also improve loadings. If a section showed trends in loss but little effect on 
loading, loss was favored in setting the leading edge. The resulting leading edge angles are 
shown in Figure 4.3.4-1 as inc_l. Relative to the baseline (Case 5), the leading edge is 
overcambered 3.5 degrees at the hub and shroud and uncambered 3.5 degrees from 55% to 75% 
span. The inc_l losses and exit swirl are shown in Figures 4. 3.4-3 and 4.3.4-4. 

The high-loss bump around 90% span was largely insensitive to incidence. Particle traces 
completed in the DAWES solution (Figure 4.3.4-5) showed that the high loss was coming from 
flow from the shroud boundary layer, the vane suction surface boundary layer, and the suction 
surface-shroud intersection. This flow was migrating down in span probably due to the vane lean 
and was unaffected by incidence at 90% span. This spanwise flow would be difficult to 
eliminate and may not be a problem relative to other stators. A similar loss source is probably 
present in all stators but isn’t so distinctive since it would stay closer to the shroud for unleaned 
vanes. 
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QHSF Stator Incidence FF 
inc_1 Exit Pt/Ptin and Particle Traces 


Figure 4.3.4-5. Stator Loss Bump 


4.3.5 Trailing Edge Angle Setting 

The vane trailing edge angle was adjusted to turn the flow as close to axial as practical. The 
trailing edge swirl angles from inc_l were used to adjust the inc_l vane trailing edge angle, 
DAWES was used to estimate the new trailing edge swirl, and within a couple of iterations the 
exit swirl was acceptably close to axial. Figure 4. 3. 5-1 shows the exit swirl for inc_l and the 
exit swirl with the reset vane trailing edge (dev_2). 


NAS A/CR— 2003-2 1 2369 


62 


QHSF Stator TE Angle Setting 
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Figure 4.3.5-1. QHSF Stator Trailing Edge Angle Setting Swirl 


The swirl is close to zero in the inner 30% span. Zero swirl in this region is important since the 
flow goes through a decreasing-radius transition duct to the core. Close control of the fan stator 
exit swirl in the core flow region is desirable to avoid mis-matching the compressor. The swirl 
farther out in span is considered acceptably close to zero. The swirl peak around 90% span is 
due to the shroud/vane intersection flow discussed previously and does not represent the actual 
stator flow turning in this region. 

Figures 4. 3. 5-2 and 4.3.5-3 show the vane trailing edge angles and the deviation. Figure 4. 3. 5-4 
shows the stator loss before and after the trailing edge angle was modified. Figure 4.3.5-4 shows 
that the loss changes caused by resetting the trailing edge angle were minor. 
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Figure 4.3.5-2. QHSF Stator Trailing Edge 



Figure 4.3.5-3. QHSF Stator Trailing Edge Angle Setting Deviation 
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Figure 4.3.5-4. QHSF Stator Trailing Edge Angle Setting Loss 


4.3.6 Composite Vane Thickness 

The vane thicknesses were adjusted in the direction normal to the suction and pressure surfaces 
to allow for composite construction. The thickness adjustment was defined by the cosine of the 
local lean angle so it was necessary to consider it last since the lean angles depend on stacking 
and on vane section stagger. 

Figures 4.3. 6-1 and 4. 3. 6-2 show stator loss and exit swirl before (dev_2) and after the thickness 
change. Meeting the composite vane thickness criteria had only minor effects on loss and exit 
swirl. Examination of the vane section loadings showed little difference. This result was 
unexpected considering how much the tangential thicknesses were increased. This may be 
because the flow was not constrained to follow 2-dimensional streamsurfaces so the effect of the 
tangential thickness increase was reduced. The “effective” vane section thicknesses are probably 
somewhere between the tangential and the normal thicknesses. Figure 4.3. 6-3 shows normal 
thickness distribution of the vane before (DEV_2) and after (THK_2d) thickness changes along 
with the recommended thickness distribution. 
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Figure 4.3.6-1. QHSF Stator Composite Thickness Loss 



Figure 4.3.6-2. QHSF Stator Composite Thickness Exit Swirl 
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Fig. 4.3.6-3. Vane Normal Thickness Distribution for QHSF Stator Vane 
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4.4 Disk / Attachment 


4.4.1 Initial Design Phase 

Since both the AE/QHSF and baseline fans will be tested in the rig, it is desirable to have the 
same fan disk for both. The primary intent during this phase was to determine the feasibility of 
using the baseline disk for QHSF with minor modifications which include changes required to 
match the disk dimensions with the NASA rig dimensions. 

The following points summarize the initial design parameter inputs for the disk design: 

TFE engine size disk scaled down to NASA rig size (Scale Factor = 22/30.698) 

• a rectangular cross section (no under cuts) 

Disk width and bore radius fixed to match the NASA rig dimensions 
(width = 3.4"; bore radius = 2.2") 

• titanium blades and attachment 

• blade airfoil load based on OPTM blade 

• blade attachment load assumed from baseline rotor configuration 


4.4.2 Cyclic-Symmetric Model 

A 2-D wedge section normal to the broach line was built using the baseline disk and dovetail 
geometry data. The bore diameter of the disk was trimmed to accommodate a 0.15" thick torque 
sleeve. The 2-D geometry of the disk/dovetail was imported to HyperMesh in IGES format for 
the finite element (FE) meshing. The FE model from HyperMesh is output in ANSYS prep7 
format for stress analysis. 

Boundary Conditions: The nodes on the two edges of the wedge were positioned at identical 
radial locations, and the corresponding nodes on two edges were coupled both in radial and 
tangential directions to enforce cyclic symmetry in the disk. The blade and disk were coupled by 
constraining blade and disk nodes at the contact plane to have same displacement in the direction 
normal to the contact surface. The blade attachment and airfoil CF loads were applied as shown 
in Figure 4.4. 2-1. The attachment geometry acts as a load transfer mechanism and was assumed 
to have zero density such that the attachment CF loads were not duplicated. The effects of slot 
broach were simulated by including a side load at the contact plane (see Figure 4 .4.2-2). Finally, 
the nodes at the bore were constrained to have zero displacement in the tangential direction. 
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4.4.3 Selection of Disk Material 


Two different materials — titanium and C250 steel — were considered for the disk during test 
runs. Analysis was performed at three points: 

• Mechanical Design Point: 100% physical speed (T t i n = 65 F) scaled to rig, 15444 RPM 

• Max. Test Speed: 1 10% Design Speed, 16988 RPM 

• Max. Operating Speed: 110% Max. Test Speed, based on NASA specification = 121% 
Design Speed, 18687 RPM 


Table 4.4.3-1. Material Properties for C250 Steel and Titanium at 75 F 


Material Type 

Iron 

Titanium 

Material Name 

Maraging Steel C250 

TI-6AL-4V 

Process/F orming 

Wrought 

Wrought 

Heat Treatment 

Double ST+Aged 

Annealed 

Yield Strength (Ksi) 

238.3 

121.5 

Ultimate Strength (Ksi) 

247.2 

129.1 

Density (lb/in 3 ) 

0.289 

0.161 


The results from 2-D cyclic symmetric analyses are presented in Table 4.4.3-2. With the 
titanium disk, the maximum radial stress was within the safe limits at the design speed. 
However, at maximum test speed, the maximum radial stress reached unacceptable limits (> 100 
Ksi) as did the burst margin. The C250 steel disk had adequate stress as well as burst margins at 
all speed conditions. Based on this analysis, it was planned to proceed with the disk design using 
C250 steel. Figure 4.4.2-2 shows the principal stress distribution in the disk/dovetail at 
maximum operating speed for C250 steel disk case. 
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Table 4.4.3-2. Maximum Stresses for Various Cases - Rig-Size Model 


Titanium Disk 


Design Speed 
(100% N phv ) 

Max. Test Speed 
(110% Np h y) 

Max Operating 
Speed (121% N phv ) 

RPM 

15444 

16988 

18687 

Max. Deflection (inch) 

0.0134 

0.0162 

0.0196 

Max. Principal (Ksi) 

87 

105 

127 

Max. Radial Stress (Ksi) 

84 

101 

122 

Max. Tangential (Ksi) 

71 

86 

104 

Av. Tangential (Ksi) 

57.8 

69.9 

84.6 

Burst Margin* 

1.42 

1.29 

1.18 

C250 Steel Disk 

Max. Deflection (inch) 

0.0105 

0.0127 

0.0154 

Max. Principal (Ksi) 

89.53 

108.36 

131.09 

Max. Radial Stress (Ksi) 

85.51 

103.47 

125.20 

Max. Tangential (Ksi) 

83.88 

101.49 

122.81 

Av. Tangential (Ksi) 

68.34 

82.70 

100.06 

Burst Margin* 

1.76 

1.603 

1.46 


* Burst Margin requirements: AE > 1.25 

NASA > 1.345 (when translated to AE criteria) 


AE and NASA have different definitions of burst margin (BM); this merits some comments. 


BM 


AE 


0.85 x F„. 


average 


> 1.25 


0.7 x F tu 

BM nasa = — > 1.5 

^ average 

where F tu is the ultimate strength of the material. 

The AE design practice recommends a material utilization factor of 0.85, whereas NASA 
recommends a value of 0.7. The AE burst factor provides a margin on rotor speed (25% over 
speed) and NASA factor provides a 50% margin on average tangential stresses. When translated 
to the over speed condition, the NASA criterion is equivalent to BM AE > 1.345 . 
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4.4.4 Attachment Stresses and Post Unzip Capability 


The margins of safety in various disk and attachment sections were obtained based on simple 
hand calculations. The sections represent the cross section that was most likely to fail under the 
given load. For example, the blade direct stresses were calculated at the blade minimum neck 
cross section. The margin of safety (MS) was defined (based on NASA guidelines) as 

Allowable Stress 

MS = TT — 7T-r~, — - 1 

Sr x Calculated Stress 

where SF = Safety Factor 

Table 4.4.4-1 provides the margins of safety for various disk and attachment sectional stresses. 
The stresses are presented at Mechanical Design point (100% N P h y ) and the margins of safety are 
presented at 100%, 110% and 121% N P h y . 

Post unzip capability represents the failure resistance of the disk post due to combined bending 
and direct stresses during blade out condition. In other words, it is the capability of the disk to 
resist unzipping or losing adjacent blades. Figure 4.4.4- 1 shows the unzip capability for the 
QHSF disk. The figure is analogous to well-known "Goodman diagram". In the figure, the y- 
axis represents a pure bending case and x-axis the pure tension case. For pure bending, failure 
occurs when bending stress exceeds 1.5 times yield strength of the material and for the pure 
tension case it occurs when the direct stress exceeds ultimate strength. For combined bending 
and direct stress case, failure occurs if the stresses are above the solid line or above the dotted 
line when safety factors are included. As seen in the figure, the QHSF disk post is safe under all 
three speed conditions. 


4.4.5 Go-Forward Design 

In summary, the preliminary design showed the baseline disk design can be used for the QHSF 
and should be manufactured with C250 steel material. It had adequate margins of safety, burst 
margin and unzip capability, and met both AE and NASA design criteria. 
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Airfoil CF load 



Side load due to 
broach angle 


Max. Principal 
Stress = 131.086 Ksi 



o 

14565 

2 9130 

43695 

58260 

72825 

87391 

101956 

116521 

131086 


Figure 4.4.2-2 Principal Stress in disk/do vetail at 121% N P h y 
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Table 4.4.4-1. Attachment/Disk Stress Safety Margins 



Stress(Ksi) 

Margin of Safety 

Percent Speed 

100% 

100% 

110% 

121% 

Contact Stress 

65.0290 

0.5630 

0.2917 

0.0676 

Blade Shear Stress 

17.8972 

2.2769 

1.7082 

1.2382 

Blade Direct Stress 

31.4828 

2.2285 

1.6682 

1.2051 

Blade Lobe Bending Stress 

14.7769 

5.8784 

4.6846 

3.6980 

Disk Shear Stress 

18.2269 

5.6909 

4.5297 

3.5700 

Disk Direct Stress 

32.5291 

5.4975 

4.3699 

3.4379 

Disk Lobe Bending Stress 

15.3263 

12.7906 

10.3972 

8.4192 

Yield Safety Factor (SF V ) = 

1.1* 


Ultimate Safety Factor (SF U ) = 

1.5* 


* Based on NASA guidelines 


Disk Post Unzip Capability Diagram 



Figure 4.4.4- 1 QHSF Disk is Designed to Withstand Unzipping During Blade- 
Out Condition 
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5.0 DETAILED DESIGN 


5.1 Rotor 

5.1.1 Geometry 

The final QHSF design configuration was chosen based on the analyses from the “go-forward” 
cases of DOE 4. Go-forward case 3 was selected as the QHSF rotor design. Final blade 
geometry span wise distributions are shown in Figure 5. 1.1-1 through 5. 1.1 -5. Tabular geometry 
data of the blade is included in Appendix I. Values in the figures represent the NASA rig rotor 
design (22 inch diameter at the rotor leading edge tip). The incidence was calculated using 
DAWES at the design point. The large incidence increase at the tip was due to the endwall 
modeling of the boundary layer at the rotor inlet. 



Figure 5.1.1-1. Blade Inlet and Exit Metal Angle 
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Figure 5.1.1-4. Blade Position of Maximum Thickness and Maximum 
Thickness / Chord Ratio 



Figure 5.1.1-5. Blade CG Stacking 
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5.1.2 Aerodynamic Results 


The design concentrated on achieving a close performance match relative to the baseline such 
that the acoustic benefits of the forward swept technology could be easily deduced from the test 
data. DAWES analyses at aerodynamic design point showed only a slight shift in speedline 
characteristics relative to the baseline. At the design efficiency goal, the flow was 0.3% lower 
than the design goal, while the pressure ratio was 1.3% higher. At the design flow goal, the 
efficiency was 0.2 points lower while the pressure ratio was 0.3% lower. However, analyses 
completed from choke to stall showed slightly higher peak efficiency than the baseline and 
approximately one point higher as the rotor was throttled up from peak efficiency to stall. 
Spanwise distributions of pressure ratio, temperature ratio, efficiency, deviation, omega-bar, and 
D-factor are shown in Figures 5. 1.2-1 through 5. 1.2-3, respectively. As previously mentioned, an 
endwall total pressure loss was being modeled at the rotor leading edge and the effect is evident 
in the figures. 



Figure 5.1.2-1. DAWES Calculated Rotor Pressure Ratio and Temperature 
Ratio at the Aerodynamic Design Point 
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Figure 5.1.2-2. DAWES Calculated Rotor Efficiency and Deviation at the 
Aerodynamic Design Point 
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Figure 5. 1.2-3. Rotor Omega-bar and D-f actor at the Aerodynamic Design 

Point 
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The increase in loss coefficient near mid-span was due to an increase in shock loss as the blade 
effective sweep approached zero degrees. Effective sweep described within is the oblique shock 
angle which the inlet flow passes coming into the passage. The oblique shock angle is a function 
of stream surface angle, blade lean angle, station lean angle (blade sweep), and relative flow 
angle. When the shock angle is considered a shock surface, both the blade leading edge and the 
impingement point on the suction surface of the adjacent blade (across the passage) must be 
accounted. Figures 5. 1.2-4 and 5. 1.2-5 show the difference in the effective sweep calculation (in 
terms of normal inlet relative Mach number) between leading edge only and a simple average 
calculation of the leading edge and the suction surface impingement point. The figures also 
show an inherent difference between the forward swept QHSF and the aft swept baseline in that 
the shock surface will always reduce the benefit of a forward swept blade, while it will always 
increase the benefit of an aft swept blade. Although the QHSF design had significant benefits 
near the tip, it was slightly worse than the baseline near the mid-span. 

Aerodynamic Design Point , NIc = 100% Altitude Op Line, N1c=90% 

Ml rel normal LE = Ml rel * sin (LE oblique shock angle) 




Figure 5.1.2-4. Effect of Blade Leading Edge Effective Sweep on Inlet Relative 
Normal Mach Number at 90 and 100% Nlc 
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Aerodynamic Design Point , NIc = 100% Altitude Op Line, N1c=90% 

Ml rel normal avg = Ml rel * sin((LE osa + SS impingement osa)/2) 




Relative Mach Number Relative Mach Number 

Figure 5.1.2-5. Effect of the Combination of Blade Leading Edge and Suction 
Surface Impingement Effective Sweep on Inlet Relative Normal 
Mach Number at 90 and 100% Nlc 


Airfoil loadings, represented as isentropic Mach number distributions along the suction and 
pressure surfaces, and the corresponding blade-to-blade Mach number contour plots are shown 
for 10, 30, 50, 70, and 90 percent spans at the aerodynamic design point in Figures 5. 1.2-6 
through 5.1.2-10. As previously mentioned, the QHSF design was chosen based on the design 
and part-speed performance results from DOE 4. In the DOE, many quality characteristics 
including airfoil loadings were examined and optimized for an overall best design. The loadings 
at 30, 50, and 90 percent spans show excessive diffusion with unnecessary re-acceleration on the 
suction surface. Additional work could be completed to improve these loadings and improve 
efficiency. 
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Figure 5.1.2-6. DAWES Calculated Isentropic Blade Loadings and Blade-to- 

Blade mach Number Contours at 10% Span at the Design Point 




Figure 5.1.2-7. DAWES Calculated Isentropic Blade Loadings and Blade-to- 

Blade mach Number Contours at 30% Span at the Design Point 
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Figure 5.1.2-8. DAWES Calculated Isentropic Blade Loadings and Blade-to- 

Blade mach Number Contours at 50% Span at the Design Point 



Figure 5.1.2-9. DAWES Calculated Isentropic Blade Loadings and Blade-to- 

Blade mach Number Contours at 70% Span at the Design Point 


NAS A/CR— 2003-2 1 2369 


83 




Figure 5.1.2-10. DAWES Calculated Isentropic Blade Loadings and Blade-to- 

Blade mach Number Contours at 90% Span at the Design Point 


Blade-to-blade contour plots are shown for 50, 70, and 90 percent spans at choke, design point, 
and near stall in Figures 5.1.2-11 through 5.1.2-13, respectively. Relative to the baseline, the 
amount of blade effective sweep significantly reduced the strength of the shock. The figures 
show that the passage shock was well contained at the tip although it moved forward as span 
decreased and was positioned upstream of the leading edge at 70% span. Attempts in DOE 4 to 
keep the shock positioned back in the passage had negative impact on rotor efficiency. 
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Figure 5.1.2-11. DAWES Calculated Blade-to-Blade Relative Mach Number 
Contours at 100% Nlc Choke for 50, 70, and 90% Spans 



Figure 5.1.2-12. DAWES Calculated Blade-to-Blade Relative Mach Number 

Contours at 100% Nlc Design Point for 50, 70, and 90% Spans 



Figure 5.1.2-13. DAWES Calculated Blade-to-Blade Relative Mach Number 
Contours at 100% Nlc Near Stall for 50, 70, and 90% Spans 


NAS A/CR— 2003-2 1 2369 


85 


A rotor-only predicted map for the QHSF is shown in Figures 5.1.2-14 through 5.1.2-16. 
DAWES analyses were completed at the points indicated. The map agrees very well with 
DAWES analyses completed with the baseline rotor. Relative to test data of the baseline rotor, 
flow and work agreed very well but the efficiency would be adjusted higher approximately 1.5 
points. 
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Rotor TR Rotor Efficiency 


0.91 


55.9% NIc 


60% 



Figure 5.1.2-15. QHSF Predicted Map (Wc vs. Eff) Based on DAWES Analyses 



Figure 5.1.2-16. QHSF Predicted Map (Wc vs. TR) Based on DAWES Analyses with 
Predicted Operating Lines 
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5.1.3 


Acoustic Results 


5.1.3.1 Design Philosophy 

The design philosophy employed in reducing multiple pure tone (MPT) noise below that of the 
baseline rotor was to use rotor forward sweep to control the strength and position of the passage 
shock emanating from the suction surface of the airfoil. It was assumed that use of forward 
sweep would delay the onset of strong shocks to higher fan speeds and allow the capture of 
shocks within the blade passages at lower fan speeds, thereby affecting buzzsaw noise as shown 
in Figure 5. 1.3-1. 


5.1.3.2 Baseline Fan 

It became apparent that the geometry and operating conditions of the baseline were such that the 
dominant shock structure consisted of a normal shock emanating from the suction side of the 
rotor, herein called a “passage” shock. At all conditions except 100% corrected fan speed (see 
Figure 5. 1.3-2), the presence of a propagating passage shock (not contained within the blade 
passage) precluded the need for a bow shock at the blade leading edge since the flow was 
subsonic at this location (see Figure 5. 1.3-3). The volume of literature gathered to date seem 
only to discuss the effects of blade geometry and stagger angle variation on the strength and 
position of a bow shock (References 2 to 13). Evidently these studies were conducted on rotors 
which operated closer to choke where the passage shock was swallowed. For this effort it was 
assumed that the geometric variations classically deemed responsible for MPT noise generation, 
due to their effect on bow shocks, in a similar way cause variations in the strength and 
propagation direction of the expelled (propagating) passage shocks. 

Using the 3-D thin-layer Navier-Stokes equation solver DAWES, the location of the passage 
shock relative to the leading edge of the following blade was predicted at a number of fan speeds 
for the baseline rotor. Figure 5. 1.3-4 shows the linear relationship between shock location and 
wheel speed. 

A correlation was developed between the predicted rotor passage shock position relative to the 
leading edge of the downstream blade and the measured buzzsaw noise in the 500 Hz, 1 kHz and 
2 kHz octave bands which contain harmonics of shaft speed below the BPF. In an attempt to 
capture noise level differences due to buzzsaw cuton, inlet sound pressures (measured from 10- 
90° from the inlet engine centerline) were summed over the 3 octave bands mentioned above and 
normalized by the level at 75.3% corrected fan speed (the speed above which MPT noise was 
cuton). Reference 14 contains measured buzzsaw noise spectra for the baseline. 
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Figure 5. 1.3-5 shows the normalized inlet noise versus passage shock position normalized by the 
rotor chord. As fan speed increases, the passage shock strength increases (see Figure 5. 1.3-6) as 
it moves closer to the leading edge of the following blade, and buzzsaw noise levels increase. 
Between 75.3% and 86.5% speed the variation of measured sound pressure level is nearly linear 
with shock location following the equation 

SPL = -45.6(x) + 13.0 

where SPL is the increase in level after MPT noise is cuton, and x is the shock location divided 
by the rotor chord. From 86.5% to 90.0% speed, however, the correlation deviates from a linear 
behavior, showing little increase in noise level with increasing rotor speed and decreasing shock 
stand-off distance. 

At 90.0% speed the shock location, predicted assuming symmetric rotor blades, is very near the 
leading edge of the downstream blade. It may be that in actuality, with asymmetric blade shapes 
and stagger angles, the leading edge of the following blade interferes with the propagation of the 
passage shock to the far-field. From Figure 5. 1.3-4 it can be inferred that near 94% speed the 
passage shock is completely swallowed and buzzsaw noise is cutoff. This is confirmed by the 
fact that in cruise (100% corrected fan speed) the baseline engine generates no buzzsaw noise 
(see Figure 5. 1.3-2). The transition between conditions of passage shock propagation and 
swallowing (as corrected speed is increased) may involve some hysterisis effects, but it appears 
to begin between 86.5% and 90.0% corrected fan speed for the baseline. 

5.1.3.3 QHSF 

As shown in Figure 5. 1.3-7, the spanwise distribution of effective sweep causes the passage 
shock location to vary with span at 100% speed. A strong shock is swallowed at 90% span while 
a weaker shock is expelled at 70% span. As fan speed is reduced to near takeoff conditions (85% 
corrected speed, Figure 5. 1.3-8), the strong tip shock dissipates leaving only a weak shock at 
70% span. 

Examination of the shock position at 2 radial locations along the rotor reveals that, at 100% 
speed and 90% span, the QHSF passage shock is severely swallowed (see Figure 5. 1.3-9). The 
streamwise shock position at 90% span is not recorded for speeds below design because the 
computed pressure gradient at that location does not suggest the presence of a shock, but merely 
a gradual pressure rise (see Figure 5. 1.3-8). At 70% span, the QHSF passage shock is expelled, 
but the strength of this shock is much less than that of the baseline at 75% speed where buzzsaw 
noise is nearly cutoff (see Figure 5.1.3-10). The shock at 90% span is relatively strong, but its 
position near the streamwise center of the rotor passage suggests buzzsaw noise cutoff. 

The main objectives have been achieved to reduce buzzsaw noise: passage shock strength has 
been greatly reduced from that of the baseline rotor and strong shocks are swallowed within the 
rotor passage. Buzzsaw noise levels are predicted to be much lower for the QHSF than for the 
baseline, if not completely eliminated. 
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Figure 5.1.3-1. Potential Buzzsaw Noise Control Accomplished Using Forward 
Sweep 



0.4 Number Contours 


Figure 5.1.3-2. Passage Shock Position for Baseline at Design Point Showing 
Swallowed Shock 
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Figure 5.1.3-3. Passage Shock Position for Baseline at Take-Off (86.5% Nlc) 
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Figure 5.1.3-4. Relationship Between Shock Location and Fan Rotational 
Speed for Baseline 
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Figure 5.1.3-5. Correlation Between Passage Shock Position and Buzzsaw 
Noise For Baseline 



75 80 85 90 

% Corrected Fan Speed 

Figure 5. 1.3-6. Predicted Strengths of Expelled Passage Shocks at the Rotor 

Corrected Speeds For Which Buzzsaw Noise Was Measured in 
the Baseline (strengths extracted from pressure distribution 
along stagnation streamline) 
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Figure 5.1.3-7. Mach Number Contours for QHSF at 100% Speed 



Figure 5.1.3-8. Mach Number Contours for QHSF at 85% Speed 
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5.1.3-9. Comparison of Shock Locations for Baseline and QHSF 
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Figure 5.1.3-10. Comparison of Shock Strengths for Baseline and QHSF 

(strengths extracted from pressure distribution along 50% pitch 
streamline) 




% Corrected Rotor Speed 


NASA/CR— 2003-212369 


94 




5.1.4 Mechanical Results 


This section presents the results of the Mechanical Analysis of the QHSF Fan Rotor Blade. 

5. 1.4.1 Design Criteria 

Mechanical design criteria were selected to satisfy both AE and NASA requirements. While the 
NASA criteria primarily addresses rig safety, the more stringent AE criteria is intended to ensure 
long-term field reliability. The QHSF blade was designed to achieve state-of-the-art noise 
reduction while maintaining commercial viability. 

• Peak steady-state stresses in the blade and attachment are at levels providing 
unlimited Low-Cycle Fatigue (LCF) life. 

• Adequate frequency margin exists to preclude resonant vibration at maximum power 
conditions. 

• Leading edge thickness and stress levels are within range of successful AE 
experience. 

• Flutter parameters involving reduced frequency and incidence are within AE 
experience. 


5. 1.4.2 Finite-Element Model Description 

A fully 3D finite element model was constructed of the entire rig scale (22-inch diameter) fan 
rotor, including the airfoil, blade attachment, and disk. The complete model is shown in Figure 
5. 1.4-1. The airfoil and blade attachment were modeled separately and connected via an AE 
multi-point constraint equation technique. Due to cyclic symmetry, the disk model consisted of a 
single, slotted wedge. Figure 5. 1.4-2 shows the detail of the blade model; the disk model will 
described in a later section. 

The QHSF fan airfoil has roughly the same geometric parameters, such as chord, hub/tip ratio, 
etc., as an existing AE fan. To reduce design costs, the platform/ blade attachment from the AE 
fan was used for the QHSF rotor. Optimization of the circumferential position of the airfoil on 
the platform was done to minimize the peak stress in the attachment. 

For computational efficiency, the disk and attachment were modeled as linear substructures. 
Blade to disk interaction at the dovetail was modeled using 3D contact surface elements. 
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Figure 5.1.4-1. Three Dimensional Finite-Element Model of the Blade, Blade Attachment, 
and Disk in Rig Scale. 
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NASA QHSF doe4 , gf 03/bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 










LD 







i — 1 




00 

o 





o 






K — 1 





o\ 



00 

ko m ko [> 

o 

c — 1 



CO 

a> 



00 

00 

H CM (J\ 

CM 



a\ 

• 

vH [> 

• O) 

2 

00 CM 

kO 00 CM 

cr\ 

CO 


in pd 

Lf) 

\ — i 

O Eh 

5 

• 

IT) 

rH • • 

kO 

• 


CN W 


kO •• 

s s 

§ 

1 

• 

• 1 

• 

kO 

N 

• fin 

01 

CM 

pq 


II 

II 

II II II 

II 

II 

II 

ii h 

>1 

; — 1 

Eh 2 

w 



Eh 




co D 

01 

pq .. 

O W 




O) 



0| 

Cd PQ 

JZ 

Woo J J 

>H 

> 

> 

> M 

pL. 



1 1 


pL, O 

w W E 

X 


ISI Q X 

>-i tsi 

5 

<C N 


' 


N 



J 


Figure 5.L4-2. Detail of the Three-Dimensional Finite-Element Model of the Blade. 
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NASA QHSF doe4 , gf 03/bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 


Material Selection and Properties 

The blade material was chosen to be Ti-6A1-4V STOA with the nominal properties summarized 
in Table 5.4. 1-1. 


Table 5.1.4-1 Ti-6AL-4V STOA Mechanical Properties 


Property 

75F 

250F 

Modulus, Mpsi 

17.5 

17.1 

Weight Density, lb/in 3 

0.161 

0.161 

0.2% Yield Strength, ksi 

137.7 

111.4 

Ultimate Strength, ksi 

146.2 

127.3 


Boundary Conditions 

A physical speed of 15444 rpm was applied to the rotor along the x-axis. The cut wedge faces 
were coupled together cylindrically to enforce cyclic symmetry. Axial (thrust) and tangential 
(torque) loads were reacted on the forward annular face representing the disk bolt flange. 

Figure 5. 1.4-3 shows metal temperatures corresponding to the aerodynamic design point as 
mapped to the airfoil model. A uniform metal temperature of 75 degrees F was applied to the 
platform/attachment and 85 degrees F was applied to the disk to simulate operating conditions. 
Design point suction and pressure-side static pressure distributions were mapped to the 
corresponding element faces on the airfoil (Figure 5. 1.4-4). 


5.1.4.3 Static Analysis Results 

As constructed, the airfoil model represented the at- speed, design point geometry (“hot” shape). 
The first phase of the static analysis is to calculate the manufactured shape (“cold”), that will, 
under design point conditions of speed, temperature, and pressure, deflect the blade into the 
desired “hot” shape. This transformation was done using an iterative procedure within ANSYS 
where the maximum error between the deflected cold geometry and the desired hot geometry is 
reduced to less than 0.001 inch. For undampered fans of this type a non-linear, large-deflection 
analysis is required. 

At the conclusion of the hot-to-cold analysis the nonlinear, steady-state stresses and deflections 
are available. Figure 5. 1.4-5 shows the maximum principal stress distribution on the pressure 
and suction sides of the airfoil. Figure 5. 1.4-6 shows the Equivalent (Von Mises) stress 
distribution. 


NAS A/CR— 2003-2 1 2369 


98 












03 







T — 1 

CO 


PO 




rH 





o 


H 


03 



CM 

00 





cr i 


s 


00 


CM 

03 

ro 

03 



00 

03 


D 


• 


CM 

LO 


O 



♦ 

^ — 1 CD 

• 

CO Eh 


03 


ro 

t — 1 

rH 

03 


p^ 

LD 

o 

o 

& *5 

LD 

CM 

T 1 

. 

• 

00 

• 


W 


o •• 

a 

a a 

o 

CM 

1 


1 

• 

CO 

N 


co 



w w 

II 

II 

II 

II 

II 

II 

II 

II 





S Ph 

a 

|Xj 


Eh 





D 

CO 

PQ .. 

o 

w s 

H 

S 


CO 




PQ 

PQ 

s 

W 03 


J w 

a 

JE] 

> 

H 



pH 


I 

<c 

o 

Ph 

W Eh 

Eh 

Eh 

£ 

Q 

X 

>H 

Csl 

p> 

Csl 








■K 

* 

* 

* 




OLnHO'vl* CO CM 00 LT) rH [> PO 
^(N^OrOLDOOHrO^OOOH^^OC^ 
(NLn<r>rOVOO^CM'X)(TiCNLnc^rNLnOO 

ro co 

• • LT)cococoi>i>i>cooocoa\a303 

LT) LD Lf) O H CN PO 'sj* LO VO |> CO 03 O ^ — I CM 

t^COCTiHHHHHHHHHHCNCNCN 


lllllllll III 


N 


i 



Figure 5.1. 4-3. Metal Temperatures Used in the QHSF Rotor Blade Static Stress Analysis. 
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NASA QHSF doe4, gf 03/bgf 03 , hot , 14194 rpm, bank temp & pres, 12/3/96 
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NASA QHSF jdoe4 , gf 03/ bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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Figure 5. 1.4-5. 


The Maximum Principal Stress Distribution on the QHSF Blade Shows 
Values Are Within Design Limits. 
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NASA QHSF doe4 , gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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NASA QHSF doe4 ; gf03/bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 



Airfoil deflections at the Design point are shown next. Figure 5. 1.4-7 shows the true-scale 
displaced shape of the airfoil (shown as element grid) along with the undeformed, cold geometry 
(shown in outline) as viewed along the tip chord. Quantitative plots of airfoil deflection are 
presented in the next four figures. Figure 5. 1.4-8 shows the vector sum (scalar) displacement of 
the airfoil. The maximum airfoil deflection was predicted to be 0.376 inch. 

Radial displacements (positive outward from center) are shown in Figure 5. 1.4-9. A maximum 
radial deflection of 0.051 inch occurred at the leading-edge. Figure 5.1.4-10 shows the airfoil 
tangential deflection (positive toward pressure-side) and Figure 5.1.4-11 shows the axial 
(positive forward) component. 

The substructure displacements at the airfoil root and the disk contact faces were used to expand 
the attachment stress and displacement solutions. Pressure and suction side attachment principal 
and equivalent stress distributions are shown in Figures 5.1.4-12 & 5.1.4-13 respectively. Figure 
5.1.4-14 shows the displacement contours for the attachment and platform. 

Inspection of the principal stress plot (Figure 5.1.4-12) and the displacement plot (Figure 5.1.4- 
14) showed a significant amount of bending in the attachment due to the large tangential sweep 
of the blade. Stress levels, however, were well within acceptable range, and the effect of the 
bending displacement was accounted for in the cold blade geometry. 


5.1.4.4 Modal Analysis Results 

Modal analysis was performed on the QHSF blade using the same model as was used for the 
static analysis (Figure 5. 1.4-2). The effect of the disk was ignored based on the relative rigidity of 
the disk combined with past experience with similar designs. 

All nodes on the contact faces of the blade dovetail were fixed in the three translational degrees- 
of -freedom (DOF). To define the Campbell diagram two analyses were performed. The first was 
done assuming uniform room-temperature, with no pressure loading, and no rotational speed. 
The second analysis assumed ADP metal temperatures (Figure 5. 1.4-3), static pressures (Figure 
5. 1.4-4), and a rotational speed of 15444 rpm. Natural frequencies for the first five modes are 
tabulated in Table 5. 1.4-2. Also included are the frequencies for the fixed-root condition 
(attachment effects neglected) under identical conditions. 
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Figure 5.1. 4-7. The True-Scale Displaced Shape of the Blade Is 


Compared With the Cold Geometry. 
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NASA QHSF doe4 , gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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Figure 5.1. 4-8. The Vector Sum (Scalar) Displacement of the QHSF Blade Is Shown With 
a Maximum Value of 0.376 in. 
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NASA QHSF doe4, gf03/bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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Figure 5.1. 4-9. Radial Displacement of the QHSF Blade Shows a Maximum Displacement 
of 0.051 in. 
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NASA QHSF doe4 , gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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Figure 5.1.4-10. The Tangential Deflection of the QHSF Blade Is Shown. 
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NASA QHSF doe4 , gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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NASA QHSF doe 4 , gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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Figure 5.1.4-12. 


Pressure- and Suction-Side Attachment Principal Stresses of the QHSF 
Blade Are Shown. I cm 
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NASA QHSF doe 4r gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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Figure 5.1.4-13. Pressure- and Suction-Side Attachment Equivalent Stresses of the QHSF 
Blade Are Shown. I cv) 
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NASA QHSF doe ~ gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 
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Figure 5.1.4-14. Displacement Contours of the Attachment and Platform for the QHSF 
<h Blade Are Shown. 
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NASA QHSF doe4 , gf03/bgf 03 , cold, 15444 rpm, bank temp & pres, 12/23/96 




Table 5.1.4-2 QHSF Vibration Results - Natural Frequencies 



Mode 

Conditions 

1 

2 

3 

4 

5 

Room temp., 0 rpm 

172 

520 

879 

1136 

1521 

ADP temp, 15444 rpm 

349 

704 

984 

1485 

1568 

Fixed-Root, 0 rpm 

177 

599 

879 

1333 

1556 

Fixed Root, 15444 rpm 

373 

791 

1020 

1600 

1682 


The resulting Campbell diagrams are shown in Figure 5.1.4-15 along with the first six excitation 
orders. Vertical lines highlight the design speed of 15444 rpm as well as the maximum test speed 
of 110 percent. Design speed frequency margin with respect to the nearest excitation order for the 
first five modes is shown in Table 5. 1.4-3. 


Table 5.1.4-3 QHSF Frequency Margin Summary 



Mode/Engine Order 

Conditions 

1/2E 2/3E 3/4E 4/6E 5/6E 

ADP temp, 15444 rpm 

58% 15% 5% 6% 2% 


As the Campbell diagrams illustrate, a significant frequency shift was experienced for modes 2 
and 4 after adding the attachment to the model. The shift in mode 2 frequency of 1 1 percent was 
especially relevant as it moved the 3E crossing to within the operating range. 

Figures 5.1.4-16 through 5.1.4-20 show contour plots of the first 5 modes at design point 
conditions with mode shapes normalized to unit generalized mass. Normalized vibratory stress 
(maximum absolute value principal stress) contours for each mode are shown in Figures 5.1.4-21 
through 5.1.4-25. 
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Figure 5.1.4-15. Cambell Diagram for the QHSF Rotor Blade Shows the Mode 3 Crossing 
in the Fan Operating Range. 
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Percent Speed (100=15444 rpm ADP) 
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Figure 5.1.4-16. First Vibration Mode of the QHSF Rotor Blade. 
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NASA QHSF, gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 2/7/97 
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*-* Figure 5.1.4-17. Second Vibration Mode of the QHSF Rotor Blade. 
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NASA QHSF, gf 03 /bgf 03 , cold , 15444 rpm, bank temp & pres, 2/7/97 
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Figure 5.1.4-18. Third Vibration Mode of the QHSF Rotor Blade. 


NAS A/CR— 2003-2 1 2369 


116 


NASA QHSF, gf03/bgf 03 y cold, 15444 rpm, bank temp & pres, 2/7/97 
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Figure 5.1.4-19. Fourth Vibration Mode of the QHSF Rotor Blade. 
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NASA QHSF, gf03/bgf 03 , cold, 15444 rpm, bank temp & pres, 2/7/97 
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Figure 5,1.4-20. Fifth Vibration Mode of the QHSF Rotor Blade. 
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NASA QHSF, gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 2/7/97 






o 




o 




00 




\ — 1 

M 



rH 

> 


in 


o 


O 



Eh 




5 


LD 


+ 





D 



(N 





w 

ro 

a> 






• 


rH 

LO 


in 

• 

rH 


• 

o 



<T\ 


• 

O 

o 

LD 

LD 


o 

o 

CO 





00 

O 

o 

r — 1 


ro 

• • 

z 


t — i 

rH 

ro 


00 

• 

i — i 

• 

CO 

rH 




II 

II 

II 


II 

II 

II 

II 



ro 

Eh 

c 

PQ 


a 

Eh 



g 

CO 

PQ 

• • 

o 

Q 

w 

PQ 

w 

S 

•* 

S 

£*1 

2 

W 

CN 

J 

O 

Eh 

tD 

os 

M 

Jg 

9 


§ 

5j 


rH 

Oh 

S 

CO 

CO 

pq 

CO 

Q 

CO 

CO 

CO 


IT) 

LD CO ^ H h rO (N CO Lf) H h ^ 

^CNroOM^Ohroo^roo^o^ 
inhroo^roo^oroo^roo^fo 

O vo 

O Toovomovoroovoroovo^o 

•VOHCNCNrO^^LnVDVOhOOCOCTiH 


111111111111111 




Figure 5.1.4-21. Stress Countours for the First Vibration Mode of the QHSF Rotor Blade. 
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NASA QHSF , gf 03 /bgf 03 , cold, 15444 rpm, bank temp & pres, 2/7/97 
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NASA QHSF, gf03/bgf 03 , cold, 15444 rpm, bank temp & pres, 2/7/97 
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NASA QHSF, gf 03 /bg£03 , cold, 15444 rpm. bank temp & pres, 2/7/97 
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NASA QHSF , gf 03 /bgf 03 , cold, 15444 rpm. bank temp & pres, 2/7/97 



5. 1.4.5 Flutter Parameters 


AE criteria for flutter is based on reduced frequency (fc/V), where f is the blade natural 
frequency, c is the chord at 75 % span, and V is the design-point relative velocity at 75 % span. 
The twist-to-flex ratio is used to characterize the mode as being flexure-dominate (small ratio) or 
torsion-dominant (large ratio). 

The reduced frequency flutter parameter calculations are shown in Table 5. 1.4-4. 


Table 5.1.4-4 QHSF Flutter Parameter Summary 


Parameter 

Mode 

f (Hz) 

c (in) 

V (ft/sec) 

fc/V 

Twist/Flex 

75% 

95% 

1 

349 

4.31 

1435 

0.6 

0.45 

0.37 

2 

704 

4.31 

1435 

1.1 

4.8 

1.1 

3 

984 

4.31 

1435 

1.6 

1.3 

1.8 

4 

1485 

4.31 

1435 

2.3 

1.2 

0.62 

5 

1568 

4.31 

1435 

2.5 

0.82 

0.42 


As can be seen in the above table, mode 1 was predominantly flexure, and the reduced frequency 
of 0.6 exceeded AE’s minimum design criteria for bending mode. Modes 2 and 3 both displayed 
significant torsional activity hence AE’s minimum design criteria for torsion mode was applied. 
Mode 2 appeared to be moderately aggressive, due in part to the frequency shift associated with 
the attachment. 

Dynamic response measurements, either with strain-gages or by other means, are highly 
recommended to ensure flutter free operation. 


5.1.4.6 Bird Ingestion 

A preliminary Foreign-Object-Damage (FOD) assessment performed on the QHSF blade 
indicated leading-edge thickness parameters to be well within successful AE experience. This 
analysis is comparative in nature and considers only spanwise geometric characteristics, LE 
thickness distribution, blade count, metal angles, bird weight, as well as other parameters. The 
computed damage tolerance factors are compared with a database consisting of very successful, 
marginal, and poor designs. This preliminary analysis however, can not account for the high 
degree of forward sweep and tangential lean built into the QHSF, which will be accounted for 
using NOSAPM. 
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A more accurate prediction of ingestion damage was obtained using the NOSAPM code 
developed for the Air Force Wright Aeronautical Laboratories (Reference 15). NOSAPM can 
perform an inelastic, transient, impact response analysis on the rotating blade, and predict with 
good success, the resulting permanent deformations of the airfoil. A key characteristic of the 
code is a loading model which allows an interactive determination of the pressure distribution on 
the impacted blade based on the blade’s instantaneous, deflected shape. With simple user input 
defining the geometry and initial conditions, i.e., rotor rpm, bird weight and trajectory, etc., the 
program solves for the blade transient response. 

Figure 5.1.4-26 shows the impact model in detail. Since the impact analysis is nonlinear, the 
model was constructed using the “cold” geometry as described above. The shaded region 
represents the predefined set of impacted elements centered upon the designated impact radius. 

As certification requirements vary with the inlet area of the fan, this analysis was performed at 
baseline fan scale. The bird weight, impact radius, as well as the rotor rotational speed were 
chosen to match conditions defined for a FAA certification test of the baseline fan. These 
conditions were defined as follows: 

Rotor Speed 88.8% 

Bird Weight 1.5 pounds 
Aircraft Forward Speed 138 knots 
Impact Radius 56% span 

Transient response of the blade due to the bird impact is illustrated in Figure 5.1.4-27 with time 
histories of the tip leading and trailing edges. The maximum displacement of approximately 3.1 
inches occurs at t=0.011 seconds at the tip trailing edge. The maximum displacement is 
overlayed on the undeformed blade in Figure 5.1.4-28. The undeformed blade shown includes 
the effects of rotational speed. 

At the conclusion of the transient (t=0.053 seconds) the permanent deformations are readily 
visible. Figure 5.1.4-29 shows the deformation contours overlayed on the undeformed blade. 
Damage appears to be localized to the leading edge tip where the maximum deformation is 0.36 
inches. The remainder of the blade appears to have been unaffected by the impact. Figure 5.1.4- 
30 presents another view of the deformation superimposed on the undeformed geometry. 

An identical analysis performed on the baseline blade revealed contrasting results. Instead of 
localized deformation like the QHSF, the baseline blade exhibited an overall gross, albeit very 
small, torsional deformation (restagger). 
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Figure 5.1.4-26. Impact Model for Assessing the Bird Strike Capability of the QHSF 
Rotor. 
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Figur e 5.1.4 -27. Transient Response of the QHSF Rotor Blade to Bird Impact. 
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Figure 5.1.4-28. Maximum Displacement of the QHSF Rotor Blade Due to Bird Impact. 
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Figure 5.1.4-29. Permanent Deformation of the QHSF Rotor Blade After Bird Strike 
Compared to the Undeformed Blade. 
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Resultant Steady-State 
Permanent Deformation 


Figure 5.1.4-30. Tip View of Permanent Deformation of the QHSF Rotor Blade After Bird 
Strike Compared With the Undeformed Blade. 
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5.2 Stator 


The analytical tools used in the stator aerodynamic design and analysis were the AXCAPS 
streamline curvature code and the DAWES 3D viscous CFD code. AXCAPS was used to 
perform a 2D streamline solution, calculating velocity triangles based on user-specified blade 
row performance profiles and generating airfoil section coordinates for use in mechanical and 
CFD analysis. These airfoil sections are generated according to user specifications of chord, 
leading and trailing edge angles (or incidence and deviation), mean camber line angles, stacking, 
and thickness. In this design, AXCAPS was used mainly as a geometry generator since the 
extreme vane stacking was felt to make some of the 2D calculations questionable. An example is 
the disagreement in leading edge swirl angle seen between AXCAPS and DAWES in the vane 
stacking FF 2. 

The rotor used in the stator design AXCAPS model changed as the stator design proceeded. The 
stator stacking DOE was done using the “optc” rotor which was the most current available at the 
time. The rest of the design was done using the final “gf03” rotor. The rotor exit velocity 
triangles which set the stator inlet velocity triangles were taken from the rotor DAWES solution 
at the design point. Some profile smoothing and adjustment was needed near the endwalls since 
the AXCAPS streamline curvature code cannot impose a no-slip condition at flowpath and airfoil 
surfaces. The flow used in the stator design AXCAPS cases was the flow predicted by the rotor 
DAWES model. 

The actual stator performance calculations used during the stator design were performed with the 
DAWES code. The DAWES calculation grid used for this stator was as similar as possible to the 
grid used in the baseline stator DAWES model, with 41 nodes pitchwise, 71 nodes span wise, and 
121 nodes streamwise. 

The model was run through enough timesteps to reach acceptable convergence. One of the 
boundary conditions in the DAWES model is the exit static pressure which, if the inlet boundary 
conditions are held constant, sets the flow. For the design point runs in the latter part of the 
design the model exit static pressure was adjusted to match the flow predicted by DAWES in the 
stator model to the flow predicted by DAWES in the rotor model to within 0.5%. 

5.2.1 Stator Geometry 

The final vane geometry was reached when the composite vane thickness criteria were met. This 
geometry is tabulated in Appendix I. 

The vane leading and trailing edge, stagger, and camber angles are shown in Figure 5. 2.1-1. The 
stator solidity, which is the same as the baseline stator, is shown in Figure 5.2. 1-2. Vane 
maximum thickness/chord, where maximum thickness is defined as tangential thickness, is 
shown in Figure 5.2. 1-3. The increase in thickness around 90% span was needed to maintain a 
desired normal thickness with an increasing vane lean angle. The waves in the thickness were 
not completely understood, but they were necessary to meet the required normal thickness 
profile; and, similar waves were seen in the required tangential leading and trailing edge 
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thicknesses. They were at least partially due to the fact that the vane lean angle at the maximum 
thickness location and at the trailing edge was slightly wavy. 

The maximum thickness location is shown in Figure 5.2. 1-4. The maximum thickness was 
located at 60% of the mean camber line length as in the baseline vane. This was set by 
mechanical, aerodynamic, and producibility considerations. The leading and trailing edge radii 
(tangential), shown in Figure 5.2. 1-5, are the radii of semicircles that form the vane section 
leading and trailing edges. They were scaled appropriately for the 22-inch rig. Edge thickness 
was approximately twice these radii. This was not exactly true for the leading edge since the 
vane leading edge was formed by a cubic generated using the semicircular leading edge as a 
starting point. The edge radii show some waviness similar to the maximum thickness. These 
were what was necessary to match the normal edge thicknesses required for composite 
producibility. 



Figure 5.2.1-1. QHSF Stator Angles 
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The vane chord, scaled for the rig, is shown in Figure 5.2. 1-6. The chord is determined from the 
solidity and vane count, which were maintained from the baseline stator. Figure 5.2. 1-7 shows 
the vane section leading edge tangential shifts applied to produce the required vane tangential 
stack. Figures 5.2. 1-8 and -9 show the chordwise distribution of mean camber line angle in 
dimensional (Fig. 5.2. 1-8) and non-dimensional (Fig. 5.2. 1-9) form. 

These plots show how the vane camber was distributed from the leading edge (0) to the trailing 
edge (1). Figure 5.2. 1-8 shows actual angles. At the leading and trailing edge, these are the hub, 
midspan, and tip angles from Figure 5. 2. 1-1. This figure provides a view of the actual curvature 
for different areas of the vane. Figure 5.2. 1-9 shows the non-dimensionalized version which is 
how these curves are usually used in design. This allows the leading and trailing edge angles (or 
incidence and deviation) to be varied without the need to change these curves. The non- 
dimensional curves were retained from the baseline stator. The curves in Figure 5.2. 1-8 were 
slightly different from the baseline stator since the leading and trailing edge vane angles were 
different. 



Figure 5.2.1-6. QHSF Stator Chord (Rig Size) 
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Figure 5.2.1-9. QHSF Stator Mean Camber Line Angles 


5.2.2 Aerodynamic Results 


Stator Streamline Curvature Velocity Triangles 

After the final vane geometry was obtained, the design point loss estimated by DAWES was used 
in the AXCAPS model along with the design point rotor pressure ratio and temperature ratio 
profiles and flow predicted by DAWES. This procedure was done to calculate performance 
information that cannot be easily extracted from DAWES results. 

Stator loss and D-factor are shown in Figures 5. 2.2-1 and 5 .2.2-2. This loss was the same as 
shown in Figure 4. 3. 6-1 except for some local modification at the endwalls necessary for 
AXCAPS convergence. The bump in D-factor around 25% span was caused by the flow splitter 
downstream of the vane trailing edge in AXCAPS. The splitter causes the streamtube area to 
increase as the flow moves radially inward and outward to get around the splitter, resulting in a 
local region of low stator exit velocity. The splitter was not modeled in DAWES. The D-factor 
peak at 90% span was caused by the bump in the loss profile from DAWES. Since this loss 
bump was due to radial migration of corner flow, the D-factor peak did not really indicate locally 
high stator loading. All D-factor values were in the acceptable range. 
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Figure 5. 2. 2-3 shows leading edge swirl from the AXCAPS calculation. Since AXCAPS 
distributes flow radially differently than DAWES, the swirl from DAWES was slightly different. 
In the AXCAPS model, the small amount of swirl at the stator exit was neglected and stator exit 
flow was assumed to be axial; therefore, Figure 5. 2.2-3 also shows the flow turning. 

Figure 5. 2. 2-4 shows the leading and trailing edge Mach number profiles from AXCAPS. These 
were slightly different than the DAWES profiles, particularly at the exit where DAWES does not 
include the effect of the splitter. These Mach numbers were within previous AE experience at 
the stator hub leading edge and at the core inlet region of the stator exit. 

Figures 5. 2.2-5, -6, and -7 show profiles of pressure ratio, temperature ratio, and efficiency at the 
stator trailing edge. These profiles combine the DAWES predictions of rotor performance and of 
stator loss and reflect the performance of the stage as a whole. 



Figure 5.2.2-3. QHSF Stator Design Point LE Swirl and Turning 
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Stator DAWES Analysis 


Figures 5. 2.2-8 through -12 show design point Mach number contours on the vane-to-vane 
surface at 10%, 30%, 50%, 70%, and 90% spans. Figures 5.2.2-13 through -17 show stator 
loadings as isentropically-calculated pressure and suction surface Mach numbers at 10%, 30%, 
50%, 70%, and 90% span. There were two things considered when evaluating these loadings: 

1. ) Judging blade row quality by looking at suction and pressure surface Mach numbers on a grid 

surface assumes that flow mostly follows that grid surface. This is not strictly the case for 
highly three-dimensional blade rows. 

2. ) The isentropic Mach number is a convenient quantity to make it easier to compare the results 

of viscous CFD codes to past experience with in viscid analyses. 

The loadings were biased towards choke near the hub which was a guideline used in setting the 
incidence, which past experience had shown to be a good strategy at the design point. 



Figure 5.2.2-8. QHSF Stator Design Point Mach Contours 10% Span 
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Figure 5.22-9. QHSF Stator Design Point Mach Contours 30% Span 

■ 1.2 50% Span Mach Contours 



Figure 5.2.2-10. QHSF Stator Design Point Mach Contours 50% Span 
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Figure 5.2.2-11. QHSF Stator Design Point Mach Contours 70% Span 



Figure 5.2.2-12. QHSF Stator Design Point Mach Contours 90% Span 
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Figure 5.2.2-13. QHSF Stator Design Point Loading 10% Span 



Figure 5.2.2-14. QHSF Stator Design Point Loading 30% Span 
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Figure 5.2.2-17. QHSF Stator Design Point Loading 90% Span 


The stator performance at several off-design conditions was evaluated using DAWES. The 
points analyzed were on the sea level operating line at 55.9%, 70%, and 80% fan corrected 
design speed and on the altitude operating line at 90%, 100%, and 105% speed. The 55.9% 
speed point is a representative approach point where most of the acoustic Y072 analysis had been 
performed. The 100% speed operating line point was slightly lower on the speed line than the 
design point so an additional run was completed. The stator inlet conditions for the stator 
DAWES model came from rotor exit velocity triangles predicted by the rotor DAWES model at 
these points. Figures 5.2.2-18 through -24 show stator profiles at these off-design conditions. 
Stator loss was well-behaved at part speed conditions. The loss bump increased between 100% 
and 105% speed and grew to affect more of the span. If efficiency at overspeed conditions (like 
max climb points) is important, this may need to be studied. The stator leading edge Mach 
numbers remained subsonic, even at 105% speed. The tip of the stator had an exceptionally 
large incidence swing throughout this speed range. The vane leading edge angle at the shroud 
was set as a compromise between design and part-speed shroud swirl. The stator exit swirl was 
adequately close to axial throughout. The importance of the locally high exit swirl 
accompanying the loss bump at 105% speed would have to be assessed considering the 
importance of overspeed performance. 
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Figure 5.2.2-18. QHSF Stator Off -Design Loss 



Figure 5.2.2-19. QHSF Stator Off-Design LE Mach 
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Figure 5.2.2-20. QHSF Stator Off-Design LE Swirl 



Figure 5.2.2-21. QHSF Stator Off-Design Incidence 
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Figure 5.2.2-22. QHSF Stator Off-Design TE Mach 



Figure 5.2.2-23. QHSF Stator Off -Design TE Swirl 
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Figure 5.2.2-24. QHSF Stator Off-Design Deviation 
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5.2.3 Acoustic Results 


Interaction of the fan rotor wake with the leading edge of the fan stator represents a significant 
source of noise at certain fan speeds. The noise associated with this interaction is most 
noticeable at approach conditions, when other noise sources are not dominant. Rotor-stator 
interaction tones are caused as the rotor wake impinges on the vane leading edge. The wakes 
emanating from the rotor are skewed and rotate circumferentially as they move downstream. 
Based on theory, the rotor wake/stator interaction noise occurs when the wake trace speed over 
the stator leading edge is supersonic. Thus, to minimize the interaction noise, the wake trace 
speed should be subsonic. 

If the shape of the rotor wake is known, then the wake trace speed may be 
controlled during the design of the stator. This can be accomplished through vane 
axial and/or tangential lean. Increased axial lean (or sweep) increases the distance 
between the stator and rotor, causing a longer time increment for the 
circumferential wake pattern to impact the leading edge tip relative to the hub. 
Increased tangential lean (pressure side down) increases the delta in slope between 
the vane leading edge and wake trace, thereby reducing the amount of 
simultaneous impingement. 

5.2.3. 1 Objective 

The objective was to use the QHSF fan rotor wake structure predicted by DAWES to help direct 
the design of the stator with the leading edge configured to minimize the acoustic effects of rotor 
wake/stator interaction, without appreciably degrading performance. The V072 acoustic 
prediction program was used to confirm the effect of the design on the interaction noise, and the 
stator design was evaluated with DAWES to determine the overall performance relative to the 
baseline stator. 

5.2.3.2 Approach 

Interaction of the fan rotor wake with the stator leading edge was assessed using a graphical 
technique. A contour plot of the rotor wake was produced on a surface formed by tangentially 
sweeping the stator leading edge curve (Figure 5.2. 3-1). By superimposing the stator leading 
edge curve on the rotor wake plot, the time required for the wake to traverse a given distance 
along the leading edge was computed, thereby obtaining the wake trace speed. 
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5.2.3.3 Design Philosophy 


The design philosophy employed in reducing interaction noise below that of the baseline was to 
use rotor forward sweep and stator aft sweep and tangential lean. The goal was to increase 
rotor/stator spacing to allow greater rotor wakes diffusion, and to increase the tangential lean of 
the rotor wakes relative to the stators thereby reducing wake trace speed and increasing phase 
variation along the stator span. 

5.2.3.4 Stator DOE 

Many stator configurations involving increased sweep and lean compared to the baseline were 
analyzed. It was found that leans of 30° and sweeps beyond the approximately 30° of the 
baseline stator did not provide significant acoustic benefit. In addition to the lack of acoustic 
benefit achieved by large leans and sweeps, DAWES results showed large flow separations near 
the hub for these cases. Consequently, more conservative leans were investigated. 

For stator DOE1 involving variations on a 15° leaned vane, interaction noise predictions showed 
that the compound bow (leaning with the rotation direction at the midspan and against rotation at 
the tip) of case 5 produced high levels while the straight and compound leans with rotation of 
cases 1, 6 and 8 produced benefits in some harmonics (Figure 5.2.3-2). The aero results showed 
that separation and high hub losses resulted from leaning the suction surface towards the endwall. 
Case 6 was selected as a go forward case. 

5.23.5 Final Configuration 

The final QHSF stage was predicted to be 3 to 5 dB quieter than the baseline in all but the 
forward-propagating 3*BPF harmonic (Figure 5. 2. 3-1). Recall that the 1*BPF tone was cut-off 
at this speed. Although the predicted wake velocity deficit was greater for the quiet fan, the wake 
width was also greater and may at least partially offset the effect of the increased deficit. 
Reduced interaction noise was likely due to increased phase variation across the stator leading 
edge as shown in the wake traces. The pressure contours were extracted from a DAWES 
solution and revealed that the 2-D kinematic analysis of V072 may predict greater wake leans 
than exist in a 3-D flow. However, variation between the Y072 and DAWES predictions 
appeared to be the same for both the baseline and the QHSF. 

The objective of reducing interaction noise using rotor forward sweep and stator aft sweep and 
tangential lean was accomplished. Large reductions in the 2*BPF harmonic will translate into 
significantly lower flyover effective perceived noise levels. 
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Figure 5.2.3-1. V072 Predictions for Final QHSF Rotor and Stator at 55.9% Speed 
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Figure 5.2.3-2. V072 Results for Stators with Moderate Leans, all with Baseline Sweep 
(30°) 
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5.2.4 Mechanical Results 


The static and vibration results for the final QHSF stator design are discussed in this section. As 
mentioned in section 4.3.3, two different retention schemes were recommended — the 
conventional single-band scheme and a two-band scheme as a backup in case flutter is observed 
in the rig. The two configurations have no aerodynamic impact; however, they have different 
mechanical behavior due to the change in the boundary conditions. 


Table 5.2.4-1 Effect of Retention Scheme on Modal Frequencies 


Mode 

Frequency (Hz) 
One-Band Two-Band 

% Change in 


Scheme 

Scheme 

Frequency 

1 

505.52 

579.66 

14.67 

2 

764.19 

740.23 

-3.14 

3 

861.4 

867.59 

0.72 

4 

997.33 

1086.7 

8.96 

5 

1296.1 

1381.7 

6.60 


Table 5. 2.4-1 shows the frequency predictions for the two retention schemes. The mode shapes 
for the first five modes are presented in Figures 5.2.4- 1 and 5.2.4-2, respectively, for one-band 
and two-band schemes. In both schemes, the first mode was a torsion mode with frequency for 
the two-band scheme about 14% higher than the one-band scheme. However, second mode 
frequency was lower by about 3%. The flutter parameter for one-band scheme was 0.85, and 
0.98 for the two-band scheme. The flutter parameter for baseline vane was 1.03 with taped 
construction and 0.96 with the braided construction. The Campbell diagram for two-band case is 
shown in Figure 5. 2.4.3. Stress analyses completed for the two schemes. Figures 5. 2.4-4 to 
5.2.4-7, showed all design criteria was achieved. 
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QHSF Composite Stator Vane; Case THK2 d ; One-Band 
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QHSF Composite Stator Vane; Case THK2 d ; One-Band 
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Figure 5.2.4-1. Vibration Modes for One-Band Scheme, Mode 1. 
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QHSF Composite Stator Vane; Case THK2d; One-Band 
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Figure 5.2.4-1. Vibration Modes for One-Band Scheme, Mode 5. 
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QHSF Composite Stator Vane; Case THK2d; One-Band 
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QHSf Composite Stator Vane; Case THK_2d; Two-Band 
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Figure 5.2.4-2. Vibration Modes for Two-Band Scheme, Mode 2. 
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QHSf Composite Stator Vane; Case THK_2d; Two-Band 
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QHSf Composite Stator Vane; Case THK_2d; Two-Band 
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QHSf Composite Stator Vane; Case THK_2d; Two-Band 
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QHSf Composite Stator Vane; Case THK_2d; Two-Band 


597 597 




o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

M3 


■sr 

00 

CM 

n 

CM 

O' 

N 

CM 

ro 

oo 



h- 


00 

r> 

o 

o 

CM 

MO 


o 

O' 

CM 

o 

00 

— 

ID 

r- 

00 

*— 

n 

If) 

r-v 

CM 

n 

O' 

CM 

VO 

h* 


— > 


•r— 

CM 

CM 

CM 

CM 

n 

n 

n 


hv 

00 

M3 



00 

<N 

ro 

CM 

O' 

r* 

CM 

n 

00 

O' 

00 


— » 

00 

n 

O 

vO 

CM 

MO 


o 

O' 

CM 

M0 

N 

^-i 

ID 

N 

oo 

»— ■ 

ro 

ID 

N 

CM 

n 

O' 

CM 

N 

00 

— 

or 

— 

— 

CM 

CM 

CM 

CN 

n 

n 

n 





CN 














O 












o o o 

^ ^ 

' I I 

M n i 



O 

O 

o 

o 

o 

o 

O 

o 

O 

ID 

o 

LD 

o 

in 

o 

ID 

ID 

O 

vO 


N 

(N 

00 

n 



n 

n 

(N 

CM 

— < 

— . 


( z i| ) Aouenbejj 

Figure 5.2.4-3. Campbell Diagram for Two-Band Scheme 
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QHSF Composite Stator Vane; Case THkI 2d; Single-Band 
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5.3 Disk 


5.3.1 2-D Analysis 

The preliminary design (PD) of the rig disk was reexamined with the optm3 airfoil loads. The 
results from the 2-D wedge model showed adequate burst margins and acceptable stresses in the 
disk. However, a 2-D axisymmetric analysis of the disk showed excessive roll at the rim with 
both baseline rotor and optm3 airfoil loads. It may be recalled that the PD of the disk was 
assumed to have rectangular cross section; i.e., no undercuts. 

The rim roll and the tangential stress distribution along the rim (slot bottom) were balanced by 
providing undercuts in the disk. Also, the disk width was reduced to 3.08" to conform with the 
blade attachment width. The bore radius was increased to 2.3" to accommodate a 0.25" thick 
torque sleeve. Figure 5.3. 1-1 shows the layout of the disk with undercuts and the torque sleeve 
mounted on the NASA balance. The tangential stress distribution in the disk hub at design speed 
is shown in Figure 5. 3. 1-2 for the two airfoil load cases. Figures 5.3. 1-3 and 5. 3. 1-4, 
respectively, show the variation of tangential stress and radial displacement at the slot bottom. 


5.3.2 3-D Analysis 

A solid model of the 3D disk was constructed in ANSYS using Boolean operators. The 2D disk 
profile was swept about the centerline by an angle of 360/22= 16.36°. A solid defined by the slot 
contour was then subtracted from the swept disk wedge, resulting in the final slotted disk sector; 
as shown in Figure 5. 3. 2-1. The finite element mesh of the final model consisted almost entirely 
of brick elements, with tetrahedra only in regions of low stress. Figures 5. 3. 2-2 and 5. 3.2-3 show 
the maximum principal stresses and the tangential stresses in the disk at design speed for QHSF 
airfoil loads. 
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Flowpath 



Figure 5.3.1-1. QHSF Disk Layout with Torque Sleeve Mounted on NASA 
Balance 
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731-60 Blade; Burst Margin = 1.35 



Optm3 Blade; Burst Margin = 1 .37 
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Figure 5.3. 1-2. Tangential Stresses in Disk Hub from 2-D Axisymmetric 
Analysis 
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Figure 5.3.1-3. Rig Disk Optimized to Perform Well With Both Blades 
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Slot Bottom Radial Deflection, inches 



Figure 5.3.1-4. Rim Rolling Well Balanced for Both Blades 
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Figure 5.3.2-2. Maximum Principal Stresses in the Rig-Size Disk at Mechanical 
Design Speed (15444 RPM) with QHSF Airfoil Loads 
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Super Element Run for 3D Disk, Rig Scale, CCW Rotation FLA, 15444 rpm 
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Figure 5.3.2-3. Tangential Stresses in the Rig-Size Disk at Mechanical Design 
Speed (15444 RPM) with QHSF Airfoil Loads 
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Super Element Run for 3D Disk, Rig Scale, CCW Rotation FLA, 15444 RPM 


Highlights of QHSF Disk are summarized below: 

• Disk Material is C250 steel 

• Disk/dovetail profile is a derivative of the baseline disk design 

• Dimensions modified to fit into the NASA balance and to accommodate .25" thick torque 
sleeve 

• Meets both NASA and AE design criteria 

• Disk has adequate burst margins, unzip capability and factors of safety 

• Optimized to perform well with both QHSF and baseline rotors 

• Well balanced to minimize rim roll by providing undercuts 

• 3-D stress analysis performed to qualify design 
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6.0 CONCLUSIONS 


The QHSF has been designed to meet aggressive acoustic goals while maintaining performance 
and mechanical qualities equal to or better than the baseline fan. The design was achieved 
through innovative use of advanced design techniques and close cooperation between the 
different disciplines on the team. Validation of this technology in the upcoming rig test will 
allow it to be used to meet the stringent acoustic requirements of the next generation of turbofan 
engines. 
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APPENDIX I 


COMPUTER LISTING OF AEROTHERMODYNAMIC PARAMETERS 
FOR QHSF ROTOR, STATOR, AND FLOWPATH DESIGN 


This appendix contains a computer listing of the aerothermodynamic output data for the QHSF 
fan stage design. All units are in the English system and headings on the columns are self- 
explanatory. Rotor and stator geometry and loss parameters are included in Tables A-l to A-3, 
respectively. Tables A-4 to A-6 contain rotor and stator performance and vector output from 
spinner leading edge (station 7) to stator trailing edge (station 22). 
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Table Al. Summary of Quiet High Speed Fan Design 


00 

00 

r— 


r" 

o CM 

oo 

o 

r~ lo 
cr> \ — i 



CP 



CP 



o 



1 — 1 



-p 



a 


n 

P 


d> 

T5 


-p 



Q) 

o 


£ 

CO 


rd 

1 


■H 




\ 



£ 


P, 



■H 

p 


-P 

-p 



rd 


P 

a 

CP 

P 

£ 

i — l 


o 

(U 

P 

1 — 1 

T5 

O 

4-1 

O 

-P 


e 

O 

-P 


P 

a 

CO 


P 

P 

P 


[2 

■H 


< 

1 

Q) 

Q 

CM 

p 


CM 

o 

£ 


a 

O 



P 

-P 

i — 1 

4H 

C 

rd 


-H 

P 

Q) 

O 

-H 

a 

P 

P 

p 



rd 

P 


£ 

i T> 


P 

■H 

I>i 

O 

CP 

P 

4-1 

<D 

-P 

P 

Q 

(U 

(D 


£ 

P, 


o 


P 

(U 

P 

rd 

tr> 

O 

P 


-P 


p> 

cd 

P 

p 

-P 

CO 

o 

CP 

HP 

a 


a 

(U 

T5 


p 

P 

P 


rd 

c 

p 



tn 

P 

< 

■H 

O 

to 

CP 

-P 


(L> 

O 


Q 

P 


>H 

P 

< 

§ 

S 

D 


oo 


P 

P 

CO CO 


P 

P 

o o 

O 

o 

O 

o o 

H 

!S 

0 


H 



\ — 1 \ — 1 

< 




P 




CO 

CO 

IP 

r- o 

CO 

CO 

P 

o 00 

< 

P 

O 

(T> 00 

P 

P 

0 


>H 

PLA 



CQ 



[2 

P 


* — 1 

* — 1 












O 

P 


00 

00 












P 

o 


CO 

CO 












P 

o 




















1 

P 

P 


U0 

U0 











1 1 

S 

CO 


5 1 

5 1 









p 


1 1 

P 

H 


CM 

CM 







PC 


uo 


1 1 

Eh 

P 










PC 


00 


1 1 












o 


LO 


1 1 


>H 










u 


5 1 


1 p 

1 

o 


LO 

00 











1 p 

H 

p 


* — i 

00 






s 





1 ^ 

P 

p 









P 





1 p 

P 

H 


00 

U0 






P 





1 p 

P 

u 


00 

00 







CO 




1 > 












P 




1 O 












p 


LO 


1 1 

CO 

o 


r-~ 

5 1 







p 


* — i 


1 1 

CO 

H 


00 

* — 1 











1 1 

p 

Eh 


00 

00 











1 1 

p 

< 








H 




LO 


to 1 

p 

PC 


5 1 

5 1 



00 


P 

u 

is 


00 


CO 







P 

o 

00 

P 

P 

o 




< 







H 

o 

00 

P 

Oa 

p 


CM 


p 1 

p 

P 


U0 

U0 


Eh 



H 

0 

p 


^P 


P 1 

s 

CO 


\ — 1 

\ — 1 



\ — 1 

\ — 1 







PQ 1 

p 

H 


CM 

CM 

O 










l l 

Eh 

P 




H 





Q 

Q 




l l 






Eh 

P 

5 1 

00 


P 

P 




l l 







C 

00 

LO 


H 

P 




l l 

1 

P 


5 1 

(T) 

P 

P 

<p 

<p 


0 

P 




l P 

H 

o 


00 

CT> 


s 




P 

CO 


^P 


1 0 

P 

p 




S 





PC 



r~ 


l < 

P 

p 


00 

U0 

> 





P 

p 


^p 


1 EH 

P 

H 


00 

00 






O 

H 


i — i 


1 to 


u 





pq 


CO 


u 

Eh 




1 1 







D 

CM 

■nT 







1 1 




CTi 

00 


P 

<p 

CM 

1 

P 




00 

1 1 

CO 

o 


00 

\ — 1 





1 

P 

pq 

0 


CO 

1 1 

CO 

H 


00 

00 




\ — 1 

o 

< 

D 

PQ 


O' 1 

1 1 

p 

Eh 








p 

> 

P 

< 



1 1 

p 

< 


5 1 

5 1 











1 1 

p 

PC 






o 

5 1 

P 













PLj 


<p 

0 

PC 



00 



Eh 



\ — 1 



H 


\ — 1 

<] 

o 



CM 



P 

PC 

IS 

* — 1 



Eh 




Eh 

P 

P 

^P 



P 

PC 

o 



>H 


\ — i 

CM 

1 

O 

H 

P 




P 

o 

p 

00 


Eh 




1 

PC 

Eh 

PC 

i — 1 



H 

0 

p 

r-~ 


H 

P 

00 

<P 













Q 

< 


O 













H 

P 

00 

r-~ 


Oa 

P 


CM 

CM 

1 

P 

p 


<p 

<p 

P 

s 



1 

s 

CO 


\ — 1 

\ — 1 

1 1 

s 

CO 


<p 

<p 

o 


\ — 1 

\ — 1 

1 

P 

H 


CM 

CM 

1 1 

p 

H 


\ — i 

\ — 1 

to 




1 

Eh 

PC 




1 1 

Eh 

P 





pq 


5 1 

1 






1 1 







D 


00 

1 


P 




1 1 


P 





P 

00 

00 

P 

1 

u 


CO 

00 

1 p 

1 

u 


5 1 

LO 





P 

H 

p 


* — 1 

^p 

1 p 

H 

p 



00 



CM 

CM 

< 

P 

p 




1 <j 

P 

p 








P 

P 

H 


o 

[> 

1 p 

P 

H 


CO 

00 





P 

P 

u 


<p 

00 

1 p 

P 

u 


<p 

O) 



00 

r-~ 

> 






1 > 







PLj 



O 






1 o 







H 

00 

00 

1 

CO 

o 


^p 

00 

1 1 

CO 

o 


o 

CM 


Eh 



1 

CO 

H 


^p 

* — 1 

1 1 

CO 

H 


U0 

CM 





1 

P 

Eh 


00 

00 

1 1 

p 

Eh 


00 

00 

0 




1 

PC 

< 




1 1 

p 

< 




P 

P 

U0 


1 

p 

PC 


1 — 1 

1 — 1 

1 1 

p 

PC 


5 1 

5 1 

H 

< 









p 






P 

P 

5 1 

5 1 







p 






P 

S 

CM 


1 

p 

P 


CM 

CM 

O 1 

p 

P 


<p 

CT> 

D 




1 

S 

CO 


* — 1 

\ — 1 

0 1 

s 

CO 


<p 

CT> 

Eh 




1 

P 

H 


CM 

CM 

1 1 

p 

H 


\ — 1 

\ — 1 


pq 

\ — 1 

o 

1 

Eh 

PC 




1 1 

Eh 

P 





D 



1 






1 1 







P 

O 

■nT 

1 






1 1 








CO 

LO 

1 

1 

P 


LO 

00 

1 1 

1 

P 


CO 

O 





P 

H 

u 


* — 1 

^P 

1 P 

H 

u 


<p 



P 



0 

P 

p 




1 0 

P 

p 





H 

<p 

5 1 

< 

P 

p 


o 

r- 

1 < 

P 

p 


U0 

00 


Eh 

<p 

CT> 

Eh 

P 

H 


<p 

00 

1 EH 

P 

H 


<p 

(Ti 



00 

CM 

CO 


u 




1 CO 


u 








1 






1 1 






P 




1 




^p 

00 

1 1 




U0 

i> 

O 

P 



1 

CO 

o 


^p 

\ — i 

1 1 

CO 

o 


■xr 

\ — i 

Eh 

< 

00 

CT> 

1 

CO 

H 


00 

00 

1 1 

CO 

H 


00 

00 

0 

P 

<p 

r-~ 

1 

p 

Eh 




1 1 

p 

Eh 




C 

s 



1 

PC 

< 


1 — 1 

1 — 1 

1 1 

p 

< 


5 1 

5 1 

P 




1 

Oa 

P 




1 1 

p 

PC 




1 
















P 





Eh 



U0 



Eh 



r-~ 



pq 




P 

PC 

IS 

CO 



P 

PC 

IS 

uo 



D 

o 

00 


P 

PC 

o 




P 

PC 

o 




P 

U0 



£ 

o 

p 

00 



S 

o 

p 

o 







H 

u 

p 

<p 



H 

u 

p 

CM 







P 

O P 


P 

O P 


P > 


P > 


P 

P 

0 

P 

p 

!S 

O 

O 

P 

p 




P 

p 

P 

H 

0 

P 

0 

Q 

< 

p 

O 

Eh 

<0 

Q 

m 




P 

P 

p 

H 

0 

0 

0 

Q 

P 

P 

O 

D 

< 

Q 

0 




P 


Q 

< 

Eh 

H 

P 

P 

S 

0 

pq 


Q 


P 

H 

>H 

< 

P 

Eh 

P 

o 

H 

s 

0 

1 


Eh 



0 

O 


P 

H 


P 

Eh 


0 

< 


< 

P 


Q 

P 


<C 

P 


O 

O 


P 

0 

p 

0 


< 


P 

Eh 

P 

P 

P 

P 

P 

P 

> 

P 

Q 

o 



p 

0 


< 



Eh 

P 


P 

P 


P 

> 


Q 

o 


<C 


O 

Eh 

p 

s 

P 

p 

Eh 

P 

p 


Eh 

p 



<c 


O 

Eh 


s 

P 


Eh 

P 



Eh 


NASA/CR— 2003-212369 


183 


ROTOR 1 6.9 5.8 .470 .418 .589 1.35 1.872 57.41 20.4 14.8 3149. 

VANE 1 8.7 8.7 .403 .403 .509 2.42 1.879 41.39 17.2 17.2 



Table A2. Quiet High Speed Fan Rotor Blade Airfoil Geometry 
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Table A3. Quiet High Speed Fan Stator Vane Airfoil Geometry 
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Table A4. Quiet High Speed Fan Rotor Blade Performance 
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Table A6. AXCAPS Design Point Data for Quiet High Speed Fan (30 pages) 
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